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1. Introduction 

In this review we present some theoretical advances in the field of quan- 
tum magnetism in two-dimensional (2D) systems. The spin-i next neighbor 
2-dimensional Heisenbera model on Bravais lattices (square, triangular) is 
Neel ordered'^ at T = OBa Frustration, small coordination number, com- 
petition between interactions can lead to specific quantum phases without 
magnetic long-ranged order. Since a decade this subject is a highly debated 
issue in the field of quantum magnetism. It was revived by the discovery of 
high-Tc superconductivity in the doped cuprates and fueled by numerous 
experimental studies of 2D antiferromagnetic insulators.B 

Section ^ is devoted to the first academic model of quantum frustration: 
the spin-i Heisenberg model on the square lattice with first- and second- 
neighbor interactions (J1-J2 model). This model is one of the most studied 
in the field and this section is a short guide to the literature, with a special 
emphasis on the various methods used for this problem. 

Section || deals with general properties of valence-bond crystals (VBC) 
and related states, the simplest phase which is commonly realized in frus- 



''This generic kind of order, with a macroscopic sublattice magnetization is called in 
the following magnetic long-ranged order (LRO), in contrast to other ordered phases 
where the long-ranged ordered correlations concern S = scalar observables (on dimers, 
quadrumers...) 
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trated spin systems without magnetic LRO. 

In section H we present large- generalizations of the Heisenberg model. 
This approach was extensively developed by Read and Sachdev from the 
early 90's and has been the first to give an insight into the alternative 
between VBC and related phases, which have long ranged order in local 
singlet patterns (whence the name of crystals), and resonating valence- 
bond (RVB) spin-liquids (SL) which are translationally invariant phases 
where the quantum coherence is a central issue. 

Section H presents some results of quantum dimer models (QDM). These 
models are effective approaches to the quantum phases of antiferromagnets 
which are dominated by short-range valence-bonds (or singlets) . They have 
received recently some special attention and provide useful insights onto 
the phenomenology of VBC and RVB SL. 

In Section ^ we review some results concerning models with multiple- 
spin exchange (MSE) (also called ring exchange) interactions. These inter- 
actions are now recognized to be present in several physical systems and 
appear to play an important role in the stabilization of RVB liquid ground- 
states. 

The last section is devoted to the Heisenberg model on the kagome lat- 
tice (and related models). Despite of an important activity on this subject, 
the understanding of the low-energy physics of the spin-i kagome anti- 
ferromagnet remains a challenging problem and we discuss some of the 
important results and questions. 

We should warn the readers that thiajjeview is quite "inhomogeneous" 
and cannot, of course, replace textbooks D^Elo^Q While some parts deal with 
some rather recent works (QDM or MSE for instance), some others are 
devoted to older results which we think are still of importance for current 
research (beginning of the section J1-J2, large- A^). The final part devoted 
to kagome reflects our own views and some unpublished material on still 
unsettled issues. Some parts are intended to be more pedagogical and con- 
crete (QDM and beginning of large- section) while some others contain 
more qualitative discussions of the physical issues (end of the section J1-J2, 
VBC, kagome). 

2. Ji J2 model on the square lattice 

We consider the following Heisenberg model on the square lattice: 




(1) 
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where (ij) and {{ij)) denote pairs of nearest and next-nearest neighbors 
respectively. Although quite simple in appearance, this spin model realizes 
several interesting phenomena which are relevant to a large class of 2D 
frustrated quantum magnets: classical degeneracy, order by disorder, de- 
struction of some long-ranged order by quantum fluctuations, break down 
of the spin-wave expansion, opening of a spin gap and (possibly ?) spon- 
taneous translation symmetry breaking, etc. For this reason we start with 
a general overview of some important results concerning this system. We 
will focus on the properties related to frustration. A review on the non- 
frustrated model ( J2 = 0) can be found in Ref.lll 

2.1. Classical ground- state and spin-wave analysis 

It is easy to find some classical ground-state of a translation invariant 
Heisenberg model on a Bravais lattice because the energy can always be 
minimized by a planar helix 

Sr = 6*1 cos(q ■ r) + 62 sin(q • r) (2) 

providedJ;hat the pitch q minimizes the Fourier transform J(q) of the 
coupling.B In the case of the J1-J2 model one has 

J (q) = 2 Ji (cos(g^) + cos{qy)) + 2 J2 {cos{qx + Qy) + cos{qx - Qy)) (3) 

• J2 < 0.5 Ji: J(q) has a single minimum at (tt, tt). It corresponds to 
the "usual" Neel state. 

• J2 > 0.5 Ji: J(q) has two isolated minima at (0, tt) and (tt, 0). They 
correspond to ferromagnetic lines (resp. columns) arranged in an 
antiferromagnetic way. These states are sometimes called collinear 
(in real space) . From these planar helix states one can build many 
other ground-states by rotating globally all the spins of one sub- 
lattice with respect to the other. Although this costs no energy for 
classical spins at zero temperature, it is known (order from dis- 
order, see below) that configurations where both sublattices have 
their staggered magnetization collinear in spin space are selected 
by thermal or quantum fluctuations. 

• J2 = 0.5Ji: J(q) has lines of minima around the edges of the 
Brillouin zone. At this point the classical ground-state is highly 
degenerate : We can write H. = est + J2 X] ('^'i + 6*2 + 5*3 + ^4) 
where the sum runs over all square plaquettes and any state where 
each plaquette has a vanishing total spin minimizes the classical 
energy. 
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Even at the lowest order in 1/5*, zero-temperature quantum correc- 
tions to the sublattice magnetization (order parameter) diverge around 
J2 = O.5J1 (Chandra and DougotEl). Such large-5 approximation usually 
tends to overestimate the stability of magnetic phases, therefore this break- 
down around J2 ~ 0.5Ji is a strong evidence for the existence of quantum 
disordered phase(s) in this region of parameter space. 

2.2. Order by disorder (J2 > Ji/2) 

The concept of "order by disorder" was introduced in 1980 by Villain and 
co-workerst2l in the study of a frustrated Ising model on the square lattice. 
In this model the next neighbor couplings along the rows are ferromagnetic 
as well as those on the odd columns (named A in the following). The cou- 
plings on the even columns (named B) are antiferromagnetic. At T = the 
ground-state has no average magnetization and is disordered. This changes 
when thermal fluctuations are introduced: a i3-chain sandwiched between 
two A chains with parallel spins has lower excitations than a B chain be- 
tween two v4-chains with anti-parallel spins. This gives a larger Boltzmann 
weight to ferrimagnetically ordered states. Villain et al. have exactly shown 
that the system is indeed ferrimagnetic at low temperature. They were also 
able to show that site dilution (non-magnetic sites) selects the same ordered 
pattern, whence the name of "order by disorder" . 

A somewhat less drastic phenomenon has been observed in quantum 
systems. It is the selection of particular long-ranged ordered quantum states 
among a larger family of ordered solutions which are classically degenerate 
at T = 0. ^ Consider a spin system in which the molecular field created 
by the spins of one sublattice on the other cancels, which is the case when 
J2 > 0.5Ji. ShendenB showed that if fluctuations are included, the system 
will select states in which all spins are collinear to each other. This follows 
from the fact that (moderate) fluctuations at one site are orthogonal to 
the mean value of the magnetization at that site and the system can gain 
some magnetic exchange energy by making such fluctuations coplanar on 
neighboring sites, that is to making the spins collinear. Such a selection of 

''In Villain's model the system is truly disordered at T = and an ordered solution 
is entropically selected at finite temperature. In the quantum J1-J2 model above, the 
classical solutions can adopt various ordered patterns: quantum fluctuations select among 
these patterns the most ordered one, that is the situation with the highest symmetry 
and the smallest degeneracy. The ultimate efl'ect of these quantum fluctuations can be 
the destruction of the Neel order in favor of a fully quantum ground-state with 0(1) 
degeneracy. 
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order by quantum fluctuations (and dilution) was discussed by HenleyEl 
and appears also quite straightly in a spin-wave expansion.113 

This selection of the (tt, 0) or (0, tt) order spontaneously breaks a four- 
fold lattice symmetry. An Ising order parameter is thus generated. It takes 
two values depending whether the ferromagnetic correlatioiis^are locally ar- 
ranged horizontally or vertically. Chandra and co-workeral^ have studied 
this mechanism and predicted the existence of a finite temperature Ising 
phase transition independent of the subsequent develomnent of a sublattice 
magnetization. This result has been questioned recfiiitlyta and the transition 
has not been observed so far in the spin-i model.llBta It has however been 
confirmed bv some recent Monte Carlo simulations of the classical Heisen- 
berg model.EZI Very similar phenomena ar e DfeiSen t in the J1-J2 quantum 



Heisenberg m^ 
Melzi el a. 



on the triangular latticed 

I have studied a quasi 2D spin-i system which is believed 
to be a J1-J2 square lattice Heisenberg antiferromagnet. They found some 
evidence (splitting of NMRhnesl for a coUinear ((7r,0) or (0,7r)) magnetic 



IRlines 



ordering. Several estimates! indeed point to J2 > Ji in this com- 

pound. 



2.3. Non-magnetic region (J2 — 

Consider the two classical "Ising states" corresponding to the wave vectors 
(tt, tt) and (tt, 0). They can be taken as (crude) variational states for the 
Hamiltonian Eq. |l|. Their energies (per site) are E^^^t^ = —J1 + J2 and 
Ett.o — — J2- As discussed above, these states cross at J2 = ^Ji. However, 
one can also consider any first-neighbor singlet (or valence-bond) covering 
of the lattice as another variational state. In such a completely dimerized 
state the expectation value of the energy per site is Edimer = —jJi, which 
is below the two Ising states around J2 — Ji/2. Of course this very simple 
argument docs not prove anything since "dressing" these classical states 
with quantum fluctuations (spin flips in the Neel-like states or valence- 
bond motions in the dimerized wave-functions) will lower the energies of all 
these trial states and it is absolutely not clear which one may eventually 
win. Nevertheless, this shows in a simple way why non-magnetic states {i.e 
rotationally invariant or spin singlet) such as dimerized states can be a 
route to minimize the energy in a frustrated magnet.'^ 



introduced a general procedure to generate local and SU (2) symmetric Hamil- 
tonians for which any first-neighbor dimerized state is an exact ground-state. These 
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2.3.1. Series expansions 

High-order series expansions can be a powerful technique to investigate frus- 
trated quantum magnets. The general method to generate zero-temperature 
perturbatioiLfixpansions in quanturn many-body systems was described by 
Singh et alE3 and Gelfand et altB For instance, one can consider the fol- 
lowing anisotropic model: 

H{x) = 2Ji J2 [s^s; + A {s^s^ + sfsy)] 

<lj> 

+2J, Yl [s!s^ + x{s^s^ + sysy)] (4) 

T-L{\ = 0) is a classical Ising model which ground-state is known. The series 
expansion about the Ising limit amounts to compute expectation values in 
the ground-state |A) of 7i(A) in powers of A: 



(A|0|A) 
(A|A> 



= ^a„A", (5) 



(energy gaps, dispersion relations and susceptibilities can also be computed 
in the same approach) . The calculation of a„ requires the enumeration and 
evaluation of the connected clusters of size ~ n, whose number grows expo- 
nentially with n. Depending on the quantity O and on the model, orders 
from 7 to 20 can typically be obtained on present computers. The series 
is then extrapolated to A = 1 by standard Fade, Dlog Fade or integrated 
differential approximations, ^uch a series expansion about the Ising limit 
was done by Weihong et alx^ foL the first neighbor square-lattice antifer- 
romagnet. Oitmaa and Weihong22l extended the series to the J1-J2 model, 
where each a„ is now a polynomial in Ji and J2. The disappearance of Neel 
order in the Heisenberg model manifests itself by a vanishing sublattice 
magnetization as well as some singular behavior of the series for Ac < 1. 
The results indicate the absence of Neel long-ranged order in the strongly 
frustrated region 0.4 < J2/J1 < 0.6. Such an expansion can locate with a 
reasonable accuracy the phase boundary of the Neel ordered regions but 
provides no direct information on the nature of the non-magnetic phase. 

To study the model around J2 — Ji/2, several other expansions have 
been carried out. A dirner expansion about an exartly dimerized model was 
done by Gelfand et alM Gelfand,Ei Singh et alH and Kotov et al.M In 



Hamiltonians arc simply defined as sums of projectors which annihilate all dimer cover- 
ings. The Majumdar-GoshO chain is the simplest example of a "Klein model". 
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this approach Ji and J2 are set to zero everywhere except on isolated bonds 
arranged in a columnar way and all the other couplings are treated pertur- 
batively. At zeroth order the ground-state is simply a product of singlets. In 
these calculations tbfi dimerized phase remains stable in the intermediate 
region. Singh et aLo also performed a different kind of zero-temperature 
series expansion. They started from a model of isolated 4-spip plaquettes 
in order to check a prediction made by Zhitomirsky and UedaEJ that such 
plaquettes could spontaneously form around J2 ~ Ji/2 to produce a state 
which is invariant under 7r/4 lattice rotations. Although the ground-state 
energy they obtained is very close to the one obtained from the dimerized 
limit (within error bars of the extrapolation procedure) they observed an 
instability in the plaquette scenario (the triplet gap vanishes before reach- 
ing the isotropic square-lattice model) which suggests that plaquette order 
is not the issue (the analysis of the exact numerical spectra for 36 sites 
confirmed this resultEZl). 

Sushkov et a/@ (improved numerical results compared to Ref.0) com- 
puted the susceptibility xd associated with the appearance of columnar 
dimer order in the (tt, tt) Neel phase by a series expansion about the Ising 
limit. Such a susceptibility seems to diverge at J2/J1 — gd — 0.405 ± 0.04. 
On the other hand the disappearance of the magnetic LRO is observed 
(through the Neel order parameter or through the anisotropy in spin space 
of the spin-spin correlations) at J2/J1 — 5c2 — 0.39 ±0.02. This point could 
a priori be different from gd . In such a case the system would first break the 
7r/4 lattice rotation symmetry at gd, while magnetic LRO remains (gapless 
spin waves). Only at gc2 > gd the SU{2) rotation symmetry is restored and 
the magnetic excitations acquire a gap. From field theoretical arguments 
based on effective actions valid close to the critical points, Sushkov et alES 
argue that the proximity (or possible equality) of gd and gc2 is a general 
feature in frustrated magnets which originate from the coupling of triplet 
and singlet excitations. 

Sushkov et alSB computed susceptibility xp associated to plaquette or- 
der by an expansion around the dimerized limit, assuming that the sys- 
tem has columnar dimer LRO. The result shows a divergence of xp when 
J2/J1 ffc3 = 0.5 ±0.02. From these results Sushkov et al. suggested that 
the translation symmetry along the columns is broken down at gcs (giving 
rise to an eight-fold degenerate ground-state in the thermodynamic limit) 
before the (7r,0) — (0,7r) magnetically ordered phase appears at 5c4 — 0.6. 
This picture is qualitatively consistent with the spin-spin correlations com- 
puted in a 10 X 10 system with a density matrix renormalization group 
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(DMRG) algorithm.^ 

Due to the relatively short series (typically of order 7) involved and 
the uncertainties in the extrapolation procedures, such results should be 
confirmed by other methods but this succession of quantum phase transi- 
tions represents a very interesting scenario. We note that if the model has a 
fully symmetric liquid ground-state in some parameter range, it should be 
difficult to capture from series expansions about limits where some lattice 
symmetries are explicitly broken. 

2.3.2. Exact diagonalizations 

Exact diagonalizations have a priori no bias, and were used very early in this 
field@'@Ei Large-size computations and sophisticated finite size scaling 
analy^ are nevertheless needed to extract significant information. Schulz 
et alx^ performed extensive exact diagonalizations of the J1-J2 model for 
system sizes up to 36 sites. They analyzed the behavior of several quantities 
such as structure factors (Neel order parameter), ground-state energy, spin- 
wave velocities (obtained from the finite size corrections to the ground-state 
energy), spin stiffness and uniform susceptibility. Their analysis, including 
quantitative comparisons with non-linear sigma model predictions,E3 con- 
cluded to the absence of Neel long-ranged order in the strongly frustrated 
region 0.4 < J2/J1 < 0.6. There, they show enhanced columnar dimer- 
dimer correlations as well as chiral ones but the size effects were not clear 
enough to discriminate between short or long-ranged order for these order 
parameters. 

2.3.3. Quantum Monte Carlo 

Quantum Monte Carlo (QMC) methods have been extensively applied to 
the J1-J2 model in the low frustration_regime giving an highly accurate 
description of the Neel phase (Sandvikcj and Refs. therein). In the non- 
magnetic and highly frustrated regime a simple QMC approach is ineffec- 
tive due to the so-called sign problem. The fixed node approach is the first 
answer to this problem: the exact imaginary time propagator e~'^^ used to 
filter out the ground-state from a variational guess |'0g) is replaced by an ap- 
proximate propagator, which has the same nodes as {ipg)- The quality of the 
result depends on the quality of the nodal regions of \ijjg). Various schemes 
have been used to try to go beyond this limitation: stochastic recnnfjgura- 
tion (SorellcO), eventually associated to a few Lanczos iterations£3'EII An 
alternative method has been devised by du Croo de Jongh et a/.,c3 where 
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the guiding function is replaced by the result of a DMRG calculation. 
Both rnfithods have their own bias. Using the first of them, Capriotti and 
SorellaEZi concluded that for J2/J1 ~ 0.45 a Gutzwiller-projected BCS wave- 
function \p BC S) was an excellent guiding wave-function: 

\pBCS) =tV\BCS) (6) 
\BCS)=e^Y>(Y.kA^<^].^\Q) (7) 

where |0) is the fermion vacuum, cj-^cf^j^ creates a valence-bond on sites 
and II projects out states with double occupancy. The pairing ampli- 
tude /ij (often called gap function A^)*^ is optimized with a Monte Carlo 
algorithm in order to minimize the energy. Capriotti and Sorella gave con- 
vincing indications that their wave-function is quite accurate. The best 
variational energies are obtained in the frustrated region with a pairing am- 
plitude which mixes d^2_j^2 and dxy symmetries. In particular it reproduces 
the correct nodal structure of the ground-state in the frustrated region at 
least for moderate system sizes where the variational result can be checked 
against the exact result. This is a subtle and non-trivial information for sys- 
tems which do not obey the Marshall's sign rule as this frustrated model. 
They concluded from these results that the system probably had no long- 
ranged order neither in dimer-dimer correlations nor in four-spin plaquette 
correlations. On the other hand, du Croo de Jongh et al. argued in favor 
of columnar dimerized phase which also break the translation symmetry 
along the columns (plaquette-like correlations similar to those found by 
series expansionsEj). 

The comparison of the results of these different approaches shows that 
this problem remains a very challenging one. The model in the frustrated 
regime is probably never very far from a quantum critical point and in 
these conditions none of the available methods seems able to discriminate 
between a VBC with tiny gaps both in the singlet and triplet sectors, a 
critical phase with a quasi order in dimers and gapless singlet excitations, 
or a true SL with gaps in any sector of spin but no long ranged order 
in any observable. As we will explain in the following sections some other 
frustrated models are happily deeper in the strong coupling regime and 
exhibit quantum phases which are easier to characterize. 

After the Gutzwiller projection Aj; is no longer the observable gap.0 
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3. Valence-bond crystals 
3.1. Definitions 

Among the different quantum solutions to overcome frustration the VBC is 
the simplest scenario. In this phase, neighboring spins arrange themselves 
in a regular pattern of singlets: dimers,*^ quadrumers or 2n-mers S = 
plaquettes. The stability of this phase comes from the extreme stability of 
small S = clusters (recall that the energy of a singlet of two spins i is 
-3/4 to be compared to the energy of two classical (or Ising) spins which 
is only -1/4), and eventually from the fact that frustrated bonds between 
two different singlets do not contribute to the total energy. 

In a VBC phase there is no SU(2) symmetry breaking, no long-ranged 
order in spin-spin correlations, but long-ranged order in dimer-dimer or 
larger singlet units. Except at a quantum critical point, all excitations of a 
VBC are gapped. Depending on the lattice geometry, such a wave function 
can spontaneously break some lattice symmetry {spontaneous VBC) or can 
remain fully symmetric {explicit VBC). In a strict sense, the name VBC 
should be reserved for systems with a spontaneous lattice symmetry break- 
ing. However, since these two kinds of systems share many similarities we 
will discuss both in this section. 

When the Haniiltonian has some incquivalcnt bonds and an integer spin 
in the unit cell (even number of spin-i for instance) the system can take full 
advantage of the strong bonds and minimize the effects of the frustrating 
ones. In that case the explicit VBC is the "natural" strong coupling solution. 
One can build a simple Hamiltonian in which the bonds which arc not 
occupied by the singlet objects are turned off. The resulting model is a set 
of small decoupled clusters (dimers or larger plaquettes) and the ground- 
state is a trivial product of singlets. Importantly, this strong coupling limit 
has the same lattice symmetry as the original one. Going back to the original 
Hamiltonian no quantum phase transition is encountered when going from 
the trivial singlet product up to real interacting ground-state. Models with 
an half-odd-integcr spin in the unit cell cannot realize a VBC unless they 
spontaneously enlarge their unit cell. In these situations there is no unique 
elected position for the 2n-mers and a symmetry breaking must take place 
in order to form a VBC. Examples of these two kinds of VBC will be given 
below. 



''whence the name Spin Peierls phase sometimes given to the VBC phase. 
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3.2. One- dimensional and quasi one- dimensional examples 
(spin-^ systems) 

One of the simplest example of (spontaneous) VBGisohsei-YEd in the J1-J2 

model on the chain for J2/J1 > (J2/Ji)c ^ 0.24.00000 For 

0.5 the doubly degenerate ground-states are exact products of dimersrcZlla 

N/2 

\MG± >= Yl |(2n,2n± 1) > . (8) 



Here and in the following we call "dimer" a pair of spins in a singlet state, 
and note it: 

-i= + > |j, - > - - > |j, + >] . (9) 

For all J2/J1 > {J2/Ji)c, the ground-states are products of dimers, 
dressed by fluctuations of valence bonds, dimer long-ranged order persists 
in all the range of parameters. This model has gapful excitations which 
can be described as_pairs of scattering spin-i solitons separating the two 
exact ground-statea£3 (these fractionalized excitations are specific of the ID 
chain) . 

The Heisenberg chain with alternating strong and weak bonds (Spin 
Peierls instability) , has indeed a unique ground-state where the dimers are 
mainly located on the strong bonds. In that case, the spin-i excitations 
are confined by theJiDjderlying potential and the true excitations are gapful 
integer magnons (E3£3 and refs. therein). It is an explicit VBC. 

A two-leg ladder with AF rung exchangehas also a unique VBC ground- 
state and gajmad magnons as excitations .e3 On the other hand Nersesyan 
and Tsvelikr^rl have proposed an example of frustrated ladder, with a 
spontaneously dimerized ground-state, and gapful excitations. Excitations 
of this last model are identified as pairs of singlet and triplet domain walls 
connecting the two ground-states, they form a continuum. 

As can be seen from this rapid and non exhaustive enumeration, VBC 
ground-states are relatively frequent in frustrated one-dimensional spin- 
i models. All these systems are gapful but the excitations could be of 
different nature emerging as modes (associated to integer spin excitations) 
or continuum of pairs of excitations that could be fractionalized (it is then 
specific of one dimensional systems) or not. 
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3.3. Valence Bond Solids 

The VBS wavp -fiinr tinn was introduced by Affleck, Kennedy, Lieb and 
Tasaki (AKLT).Eile3 It can be constructed whenever the spin S' on a site 
is a multiple of one half the coordination number z: 25* = mod z. Let us 
consider the simplest case 2S = z. In that case the local spin 5" can be seen 
as the symmetric combination of 25' (fictitious) spin-i. Now on each bond 
of the lattice one can make a singlet between two fictitious spins- i. Such 
a product of singlets does not belong to the physical Hilbert space of the 
original spin-S' model but to a much larger space. The VBS wave-function 
is defined as the projection of the singlet-product state onto the physi- 
cal space. This projection amounts to symmetrize (for all lattice sites) the 
wave-function with respect to the fictitious spins to force them into a phys- 
ical spin-S* state. A VBS can be viewed as an explicit VBC of the fictitious 
spins. Simple Haniiltonians with short ranged and SU (2)-symmetric inter- 
actions for which the VBS is an exact ground-state can be constructed (sum 
of projector Si) . By construction the VBS wave-function is a spin singlet 
and breaks no lattice symmetry. By extension we may say that a system is 
in a VBS phase if its ground-state can be adiabatically transformed into the 
VBS wave-function without crossing a phase transition. As the VBC, mod- 
els in the VBS phase have a gap to all excitations'^ but their wave-functions 
are slightly more complex and their order parameter is non-local. The or- 
der of VBS is associated to long-ranged singlet-singlet correlations in the 
fictitious spins. Expressing such observable in terms of the pljys iffil spins 
leads to a non-local quantity called string order porameier.Ej^Eil Contrary 
to explicit VBC, VBS have fractionalized degrees of freedom at the edges 
of the system with open boundary conditions. These are simply associated 
to the unpaired fictitious spins. To our knowledge these properties have not 
been explored in quantum 2D systems. 

The spin-1 Heisenberg chain is the prototype of VBS in ID.'^ Such a 
state has a unique ground-state, a gap in the excitations and exponentially 
decreasing spin-spin and dimer-dimerE3 correlations. See the chapter by 
P. Lecheminant in this volume for more details about the VBS phase of the 
spin-1 chain. 



''This.^pay however not always be true when the coordination number of the lattice is 
large.tj In such cases the VBS wave-function is still a spin singlet but has long-ranged 
spin-spin correlations. We do not consider such cases here. 

'^In ID, some authors call "Haldane systems" all the spin-gapped phases, whatever their 
true nature: VBC or VBS. 
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A spin-| specific SU (2)-invariant model on the honeycomb latticeElii is 
another example of 2D VBS. The spin-f Heisenberg modeLon the kagome 
lattice was proposed to realize a VBS-like ground-stateEj in which sin- 
glets form on every hexagon without any spontaneous symmetry breaking 
(hexagonal singlet solid). A similar approach was carried out for the spin-1 
pyrochlore Heisenberg modeLEll In that case a lattice distortion was invoked 
to lift the degeneracy between the two singlet states of each tetrahedron. 



3.4. Two-dimensional examples of VBC 

3.4.1. Without spontaneous lattice symmetry breaking 

Two spin-i experimental examples of 2D (explicit) VBC have recently at- 
tracted attention: 
CaV4O9ii0i@li0 and SrCu2(BO3)2.li§fflSSilliS0 Both of 
them have a lattice derived from the square lattice: in CaV409 the V''+ ions 
are on a 1/5 depleted square lattice and SrCu2 (603)2 is an illustration of 
the Shastry Sutherland modelH . A lattice embedding the main couplings 
of these two physical problems is drawn in Fig. ^ Interactions are of the 
Heisenberg type.° 

VBC are obvious ground-states in the strong coupling limits of each 
problem. 

The exact ground-st ate o f the Shastry Sutherland model is built from 
singlets on the J' bondsSEZi For J' / J > 1.55 ±0.05 the model has a gap in 
the spectrum of excitations and is in a dimerized VBC phase. ^ For J' /J < 
1.15±0.05 the system is in the (tt, 7r)^eel state of the square lattice (results 
of zero-temperature series expansiorO) . The possibility of an intermediate 
phase, possibly with helical shor t-ranged correlations, has been actively 
discussed in the literatureEaSH S. Miyahara and K. Ueda have recently 
written a review of the theory of the orthogonal dimer Heisenberg spin 
model for SrCu2 (603)2 



•^Each kagome site belongs to two hexagons. Each physical spin-1 can be split into two 
spin-i, each of them being involved in the formation of a singlet on one neighboring 
hexagon. r—i 
■^Small Dzyaloshinsky-Moriya interactions have been identified in SrCin(RO->)T.P^I 
^Consider a decomposition of the Shastry-Sutherland lattice as edge- and corner-sharing 
triangles. For J' = 2 J the Hamiltonian is accordingly written as a sum of J {§1 + 32+83)'^ 
terms for each triangle (J'-bonds are shared by two triangles) and each such term is 
minimized by the dimerized state. This shows rigorously that the dimerized state is the 
ground-state as soon as J' /J ^2. 
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Fig. 1. The i-depleted lattice and the Shastry Sutherland lattice. The strong bonds 
of the Shastry Sutherland model are the bonds J' (dotted dashed lines): they can ac- 
commodate orthogonal dimers ( J = K acts as a perturbation in the real SrCu2(B03)2). 
The lattice formed by the sjtwjngest exchanges in CaV409 is obtained with J = 0. The 
phase diagram of this modolE3 contains collinear Neel phases, dimer and 4-spin plaquette 
VBC. 

The i-depleted Heisenberg square lattice model [J — 0) has been stud- 
ied as a function of the ratio of the two different couplings: bonds within a 
plaquette {K) and dimer bonds (J') between plaquettes. At isotropic cou- 
pling (J' = K) collinear long-ranged Neel order survives the depletion, (the 
decrease in the order parameter is roughly 50%O). A small unbalance in 
couplings drives the system either in a 4-spin plaquette VBC {K > J') or 
in a dimer VBC {K < J'). Both (exphcit) VBC phases have a spin gap. 
A recent generalization of these models by Lauchli et al encompasses both 
the -l-depleted Heisenberg square lattice model and the Shastry Suther- 
land modeS (see Fig. |]). Its phase diagram exhibits collinear Neel phases 
((tTjTt) or (0,7r)) separated from the VBC phases by second order phase 
transitions. Transition between the two VBC phases which have different 
symmetries occurs via a first order phase transition. s 

Excitations in these models come from the promotion of local singlets to 
triplet excitations. In 2D the ordered dimer background provides a confin- 
ing force for the spin-i excitations. Indeed, separating two unpaired spins 
(that is two spinous) creates disruption in the ordered pattern all the way 
from the first to the second. The energy cost is thus proportional to the 
length of the string of defaults and both spin-i excitations remains in fact 



recent Sp{N) study of the Shastry Sutherland modeO suggests that a spin liquid 
phase with deconfined spinous might appear in such a model. No evidence of such a 
phase has emerged from the SU{2) studies. 
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confined. Only integer spin excitations are expected. On the other hand in 
these strongly coupled models single-triplet hoppings can be difficult and 
correlated motions might be important, leading to a large zoology of ex- 
cited modes (see Ref.Ej and Refs. therein). This potential frustration of 
the triple t rn o tinn fa vors the appearance of magnetization plateaus in VBC 
models-olIaO'EBE^O Thisaspect was briefly discussed in the lecture notes 
published by the authors £3 



3.4.2. With spontaneous lattice symmetry breaking 

In the previous models the (explicit) VBC phases do not break any lat- 
tice symmetry. They can be directly related to the geometry and relative 
strength of the couplings. In more symmetric situations with frustration, 
spontaneously symmetry breaking VBC can appear as a way to overcome 
this frustration by taking full advantage of the quantum fluctuations. This 
is probably the case in the J1-J3 model on the square lattice for inter- 
mediate J3/J1 ^ 0.6,0 in the J1-J2 model on the hexagonal lattice for 
intermediate Jt /, 7i ^ ^-.^-^ a-nd in the Heisenberg model on the checker- 
board latticcBE '!£!^ 10^,102,103 |-]^g gj.g,(. (ja^ggg ^j^g ground-states are 
dressed columnar VBC of dimers. Translation and C4 (resp. C3 only) sym- 
metries of the lattice are spontaneously broken. The ground-state is 4 (resp. 
3) times degenerate. Spin-spin correlations decrease exponentially with the 
system size. All excitations are gapped. Cnntpr y to the J1-J2 model on 
the square lattice, exact diagonalizationsHaO^Lij give a rather straightfor- 
ward information on these systems where the correlation lengths are small 
enough (far enough from the critical points which limit the boundary of the 
VBC phases). 

The spin-i Heisenberg model on the checker board lattice, which 



be 



105 



101, IOC. 39 



as a planar lattice of corner sharing tetrahe- 
hf\ fj rece ived the largest attention for different 
123 The problem has classically a continuous 



102 



can also 
drons (se 
reasons J12 

local degeneracy: the Hamiltonian can be rewritten as the sum of the 
squares of the total spin of each tetrahedron, and every configuration 
with a zero spin on each tetrahedron is a ground-state. Classically this 
problem shar e s t his property with the Heisenberg model on the kagome, 
pyrochlorell2j'E2ZI, garnetE23 and pyrochlore slabll23 lattices (these lattices 



made of corner sharing "simpl exe s" with 2, 3 or 4 spins each were dubbed 
"bisimplex" lattices by HenleylH3). 

The quantum spin-^ antiferromagnet on the kagome lattice has been 
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^bin S=0 



Fig. 2. Left: the checkerboard lattice. The spins sit at the vertices shown by bullets, 
all couplings are identical, ui , U2 are the unit vectors of the Bravais lattice. Right: the 
ground-states of the Heisenberg model on the checker board lattice are translational 
symmetry breaking VBC of 4-spin plaquettes on the uncrossed squares. 



found to be quite specific with a small gap (if any) toward triplet excitations 
and an anomalous density of gapless low- lying singlets excitations (see §0). 

The quantum scenario on the checkerboard lattice is quite different. 
The ground-state is a (dressed) product of 4-spin 5 = plaquettes on 
uncrossed squares: this state breaks translational symmetry but not C4 
(the point group being defined at the center of an empty plaquette). It has 
a two fold degeneracy in the thermodynamic limit: this is easily seen in the 
symmetries and finite size scaling of the low lying levels of exact spectra.cSI 
The choice of the 4-spin 5 = states on the uncrossed squares corresponds 
to the most stable local configuration. The product oLjS* = quadrumers 
is only weakly renormalized in the exact ground-stateJ2j 

Excitations of thismodel have been studied in different approaches, 
exact diagonalizationalij, series expansionsJ122l real space renormalization 
transformation.1123 All the excitations (singlet and triplet) are gapped. The 



triplet excitations originate from the triplet excitation of an uncrossed pla- 
quette, they are weakly dispersive with a large gap. The singlet excitations 
cannot be described so simply: from exact diagonalizations data one can 
suspect that some of these excitations come from the reorganization of 
two adjacent triplet on crossed squa,res. The Contractor Renormalization 
(CORE) method of Berg et aZ.E£3 on the other hand suggests that these 
excitations are domain walls between the two translated plaquette- VBC 
ground-states. 
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3.5. Methods 

Spin wavesi^ and Schwinger bosonsS0 are simple approximations to study 
the phase diagram of a quantum frustrated magnet. But the first method 
only gives an approximate knowledge on the range of existence of the Neel 
phases, and it is rather difficult to include the effect of fluctuations be- 
yond the mean-field approximation within the Schwinger-Boson formula- 
tion. As discussed in the section devoted to large- approaches (|Q , spon- 
taneous lattice symmetry breaking are very likely to be due to topologi- 
cal effects (Berry phase of instantons) which cannot be captured by the 
mean-field state. This probably explains why no spontaneous VBC has 
ever been found (to our knowledge) in Schwinger-Boson calculations.'^ For 
unfrustrated problems, as the Heisenberg model on the ^-depleted square 
lattice,Lil QMC is considered as the method which can give benchmarks for 
other approaches. 

Although VBC are naturally described with the help of spin-i valence- 
bonds, the appearance of such states in the low-spin limit can sometimes be 
anticipated from an appropria te la rge-iS* approach. An example is give n in 



the work of Tchernyshyov et a/Ji2j (see also a previous paper by Henleylii3) 



on the checkerboard Heisenberg model. In this model, when both couplings 
(square lattice bonds and diagonal bonds) are equal, ground-state has a con- 
tinuous local degeneracy. However leadirig 1/S corrections select collinear 
states out of this huge manifoldJlij There remains an exponential number 
of such states and they do not have any magnetic order. However thev 



exhibit long-ranged bond order and a spontaneous symmetry breakingt2; 
which is analogous to the one observed numerically in the spin-^ case. 

For frustrated problems, exact diagonalizations can be useful tools in 
situations where the system is not too close from a critical point, that is 
when the correlations length is not too large. Successful applications of ex- 
act diagonalizations methods to 2D Heisenberg magnelja^realizing a VBC 
include studies of theJp^-Js model on the square latticej23 J1-J2 modsLon 
the hexagonal lattice,Lll Heisenberg model on the checkerboard latticeO In 
such situations analysis of the quasi degeneracy of the low- lying levels of the 
spectra and of their finite-size scaling gives an unbiased and direct informa- 
tions on the symmetry breakings in the thermodynamic limit O However 
the boundaries of the phases and the quantum critical points cannot be 



''It seems however that spontaneous VBC naturally arise in large-Af approaches with 
ferm ionic representation of SU{N) when 1/A'^ corrections are considered (see section 
0). 
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accurately determined with this method. The series expansions described 
in the previous s pctinn appears to be a powerful approach to determine 
phase boundaries .Ili3ilj If the spin-spin correlation length is not very short- 
as it is the case in the J1-J2 model on the square lattice for J2/J1 ~ 0.5,e3 
it is very difficult to decide from exact diagonalizations between a VBC, a 
critical phase or an RVB SL. 

Concerning excitations, exact diagonalizations give the gaps in each 
sector and provide a crude approximation of the dispersion laws of the 
first excitations. The large-scale nature of the excitations (as for examples 
domain walls excitations) can escape this method. The semi-analytical ap- 
proaches which can be used for the study of the excitations of the VBC, 
all use as a basic departure point the excitations of a local cluster of spins 
conveniently renormalized by perturbationE^ ]£i or effective renormaliza- 



tion techniqueJl£3 Contrary to exact diagonalizations these methods are 
not limite d by fin ite-size effects but the results can be biased by the depar- 
ture point .E 



3.6. Summary of the properties of VBC phases 

The generic features of VBC (whatever the dimensionality of the lattice) 



are: 



A spin gap, no SU (2) symmetry breaking and short-ranged spin- 
spin correlations. 

Long-ranged order in dimer-dimer and/or larger S = Q plaque- 
ttes. The coupling of this order to lattice distortions is probable in 
experimental realizations of spontaneous VBC. 
In spontaneous VBC phases the ground-state is degenerate. From 
the theoretical point of view the discrete symmetry of the order 
parameter of the VBC which spontaneously breaks a lattice sym^ 
metry may give birth to a finite temperature Ising-like transition.lU 
Simultaneity between this transition and a possible structural tran- 
sition is likely when the couplings o f th e spins to the lattice degrees 
of freedom (phonons) is considered. 



113 



VBC have gapped excitations, in the S = sector as well as in 
other S sectors. A wide zoology of modes is to be expected as well 
as continua associated to multi-particle excitations or scattering 
of domain walls (in the case of a spontaneous symmetry breaking 
of the ground-state). In two dimensions all these excitations have 
integer spins (the ordered back-ground inducing a confinement of 
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the spin-i excitations) 

Frustration on the square lattice or more generally on bipartite lattices 
is often overcome by VBC phases. No theorem forbids the appearance of 
VBC in triangular geometries but there is up to now no examples of such 
phases in pure spin-^ models (in Sec. || examples will be given within the 
framework of quantum dimer models). 

It has been advocated in the large- iV approaches (see section^) that, at 
least in two dimensions, coUinear spin-spin correlations generically lead to 
VBC or VBS and only non-coUinear spin-spin correlations can give birth 
to RVB SL with unconfined spin-i excitations. The present knowledge of 
SU{2) phase diagrams supports this prediction. The VBC found so far nu- 
merically in SU (2) spin models appear to be in regions of parameter space 
where the spin-spin correlations are characterized by some short-ranged 
coUinear order in the large-^* limit. The J1-J2 model on the honeycomb 
lattice has a classical incommensurate phase in the regime of high frus- 
tration and there are some evidences that in the quantum phase diagram 
thejx)llinear phase is separated from the columnar VBC phase by a RVB 
SLEZI The multiple-spin exchange (MSE) model on the triangular lattice is 



also believed to be a RVB SlJli^ and the corresponding classical grouniL 
states generically have non-coplanar spin configurations. Capriotti et aZB3 
argued that the spin-i square lattice J1-J2 model could be a RVB SL. If 
confirmed, this would be the first counter-example to the general rule ex- 
plained above (The Hcisenberg model on the pyrochlore lattice might be 
an other counter-examplellHa). 

4. Large- AT methods 

Introduced by Affleck ,lliiAflb£k and Marston, Li^ Arovas and Auerbachllii 



and Read and Sachdevti2llli3llH3 in the context of spin models, large- iV ap- 



proaches are powerful methods to investigate quantum antiferromagnets. 
When N is taken to infinity many of these models can be solved by saddle 
point methods and finite-iV corrections can be, at least in principle, ex- 
plored in a controlled way. A success of these approaches is that they can de- 
scribe the phenomenology of a large variety of phases encountered in quan- 
tum magnets : Magnetic LRO (possibly with order by disorder selection) as 
well as phases dominated by short-ranged valence-bonds VBC VBS and 



116 



119 



12o[l21 



RVB liquids. One crucial result (due to Read and Sachdev 
concerning these three later phases is that the analysis of finite- correc- 
tions to some large- A'' formulations {Sp{N) for instance, see below) provides 
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a general criterion to decide which ol these three phase should appear in a 
given model. 

This criterion is the following in 2D: if the (large- TV equivalent of 
the) "spin" S at each site matches the lattice coordination number z by 
25 = mod z a VBS is to be expected. If it is not the case (as for a 
spin-i model on the square lattice) one should look at the local spin-spin 
correlations. If they are reminiscent of a coUinear order, a VBC with spon- 
taneous translation symmetry breaking is expected whereas non-coUinear 
short-ranged correlation generically give rise to a RVB phase without any 
broken symmetry and deconfined spinon excitations. These results are of 
course based on a large- iV generalization of the original spin model and 
there is no guaranty at all that these rules should always apply to SU{2) 
models. To our knowledge they have however not been manifestly found in 
error up to now.' In the following we will present some of the important 
reasoning steps leading to this result. 

4.1. Bond variables 

The SU{2) algebra of a spin S at one site can be represented by 2 species 
of particles aj. (with a =t, [), provided that the total number of particles 
on one site is constrained to be aja| + aja^ = 25*. The raising operator S*"*" 
(resp. S~) is simply represented by ajaj^ (resp. ajai). These particles can be 
chosen to be fermions (Abrikosov fermions) or bosons (Schwinger bosons). 
These particles carry a magnetization ±i since S"^ = i(aja| — ajax). For 
this reason they are often called spinous. The Heisenberg interaction is a 
quartic interaction for these particles: 

■ S, = S^- \a\^A,, (10) 
with the bond operator Aij defined by: 

A,.j = ajj^a^i - a^^aii (11) 

Acting on the vacuum, A}^ - creates a spin singlet on the bond (ij). Physically 
A^jAij measures the number of singlets on that bond and Eq. |l^ shows 
that the antiferromagnetic Heisenberg interaction just tries to maximize 
that number. 

The idea of large- methods is to generalize the SU{2) symmetry of the 
spin— S* algebra to a larger group SU{N) (or Sp{N)) by letting the index 



'The spin-i kagome antiferromagnet might however be such an example. See section 
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cr go from 1 to iV (or 2N). The SU{N) (or Sp{N)) generalization of the 
Heisenberg model is solved by a saddle point calculation of the action. The 
N ^ oo limit is very similar to a mean-field decoupling of the four-body 
interaction of the physical SU{2) model: Aj^A.j ~ {A,j) + A^- - 



4.2. SU{N) 

The generalization of the Heisenberg model to such a symmetry group 
depends only on the choice of the irreducible representation of SU{N) ac- 
cording to which the "spin" operators transform (and not on the choice 
of fermions or bosons to implement the representation). For SU{2) this 
amounts to specify the magnitude S of the spin. For SU (N) irreducible 
representations are labeled by Young tableaux. The case of a general rect- 
angular iableau with He columns and m rows was discussed by Read and 
Sachdevll23 and Uc appears to continue to play a role similar to 25" in the 
large- -/V phase diagrams In this review we will focus on a slightly different 
large- TV generalization of the SU{2) model which is both able to deal with 
frustration and magnetic states. 



4.3. Sp{N) 

To perform a large- ex tension of frustrated Heisenberg models one has 
to use either fermionaHj or bosons with an Sp{N) symmetry. The latter 
seems to produce phase diagrams that closely resemble the SU{2) prob- 
lems and we will focu s on this representation which was introduced by 
Read aniLSac- 



papers 



lis 



nq . 



11*120 



idpv l^^m Thp presentation below is largely inspired from their 



We now have 2N flavors of bosons at each site: hi^a with a = 1..2N and 



j Taking the limit N ^ oo with m fixed of order 1 and Uc ^ N is most conveniently 
done with bosons fej,p where a = I - ■ ■ N is the SU{N) index, while p = 1 ■ ■ ■ m label 
the different "colors" . There are therefore Nm kinds of bosons. On the other hand it is 
convenient to use Uc "colors" of fermions (still with an SU{N) index) to deal with the 
case Uc ~ 0(1) and m N . Bosonicj'epresentations with ^ N are appropriate to 
describe magnetically ordered phasest£j but cannot be used for frustrated models (such a 
representation is not self-conjugate j i . Qn t he other hand fermionic representations, such 



as the m = N/2 and = 1 one,tijt£f| can be used on any lattice l put || th i j^r (| o not 
display magnetically ordered phases and tend to favor dimerized states 
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we define an iS'p(A^)-invariant bond operator: 

(Jm' = 1..2N 

where the 2N x 2N antisymmetric tensor J is block diagonal 



(12) 



J 



1 

-1 



1 

-1 



(13) 



and generahzes the SU{2) antisymmetric tensor {SU{2) is identical to 
Sp{l)). Up to some constant the Sp{N) Hamiltonian is 



with the constraints 



27V 



^J'^ = ^'^ 



(14) 



(15) 



ric — 2S in the SU{2) case and Uc/N — k will be kept constant when 
taking the large- limit. The partition function can be represented by an 
imaginary time functional integral: 



Z = Jv[X,,bf,bl^]e,^p(^~j\Lo+n)dT^ 
Lq = ^6|°"(5r + iXi)bia - iNn'^Xi 



(16) 



(17) 



The Xi (r) are Lagrange multipliers that enforce the constraint Eq. |5| at 
every site. Bond degrees of freedom Qij are introduced in order to decouple 
the bosons (Hubbard-Stratonovitch). The partition function is now 



Z 



V[Q.u,Q,,,X,,b'^,bl^]exp 



{Lo + Li)dT 



with 



J a 



-(4Q,,+h.c) 



(18) 



(19) 
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SO that a Gaussian integration of the Qij gives back Eq. ^ The bond 
variable are Sp{N) invariant and they can take non-zero expectation values 
at a mean-field level without breaking the Sp{N) symmetry. As we explain 
below, they are however not gauge-invariant. 

4.3.1. Gauge invariance 

An important property of Eq. |l^ is the U{1) gauge invariance associated 
to the following transformations: 

b,, ^ b,„e"^^ (20) 
Q^^ ^ Q^^eK*-+M (21) 

A. ^ A, - ^rcb^ (22) 

where 4'i{T) are arbitrary site- and imaginary-time-dependent angles. This 
gauge invariance comes from the the conservation of the local boson number 
and reflects the fact the magnitude of the spin is constant at each site. 
If we focus on the phase degrees of freedom of the bond variables, the 
Eq. |l^ describes a system of charge- 1 bosons coupled to a U{1) lattice 



gauge theory££3'E£3 These gauge degrees of freedom play a crucial role in 
the analysis of the fluctuations about mean-field solutions. 

Effective action for the bond variables — The boson degrees of freedom 
can be integrated out to give an effective action for the bond variables: 

Z = y"l?[Qy,0.j,A,]exp(-5°*f) (23) 



TVTrlogG (24) 



where is the quadratic form which couples the bosons in Eq. |T8| (prop- 
agator). It depends on the bond variables and on A^. We may write formally 
— dr — iX — Q. The term iVTrlogG is the free energy of the bosons 
in presence of the bond fields. By construction the action S'^^ is gauge- 
invariant with respect to the transformations of Eqs. |20|-|2^. So far this is 
an exact formulation of the original model for arbitrary N. 

4.3.2. Mean-field (N ^ oo limit) 

Since N factorizes (no flavor index is left in Eq. ^), Z is dominated by 
the saddle point of S'^^ when N is large. For simple models such as the 
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first-neighbor antiferromagnet on the cubic lattice (any space dimension), 
the saddle point can be determined analytically. The TV = oo limit is al- 



most equivalent to the Schwinger-boson mean- field theoryH3'E!l This can 
otherwise be done numerically.'^ In this large- TV limit, two kinds of mean- 
field solutions can appear. For large enough k the bosons condense at some 
wave-vector, the spectrum of the mean-field Hamiltonian is gapless. This 
corresponds physically to Neel long-range order. On the other hand, for 
smaller k (smaller "spin") the mean field Hamiltonian is gapped and the 
ground-state preserves the Sp{N) symmetry. Fluctuations around the sad- 
dle point are not expected to change drastically the Neel ordered phases 
but they play an important role in the physics of Sp{N) symmetric phases. 
The following is a brief discussion of the effects of fluctuations in these 
non-magnetic phases. 

4.3.3. Fluctuations about the mean-field solution 

At the mean-field level described above some Qij acquire a (static in all 
known cases) non-zero expectation value: (Qij) — Qij. For this reason such 
a state spontaneously breaks the local gauge invariance of Eqs. and 



22[ However this does not mean that the gauge degrees are all gapped 
and do not play any role at low energy. In fact, as remarked by Read and 
Sachdev, depending on the geometry of the lattice defined by the bonds 
where Qij ^ 0, some long- wavelength gapless gauge excitations survive and 
the associated fluctuations must be taken into account. More precisely, the 
fluctuations of the bond variables about the saddle point are decomposed 
into an amplitude and a phase 

Qn = + 9y) exp(i%) (25) 

and we expand S'ofr with these new variables. Two cases must then be 
considered: 

i) The lattice made of the sites connected by non-zero Qij bonds 
is bipartite. This is automatically the case if the original lattice 



''To our knowledge, all the saddle points considered so far are static (expectation values 
of the Qij are time-independent) and the corresponding (Qij) could all be made real 
with an appropriate gauge transformation. There is no chiral order and the time-reversal 
symmetry is unbroken. The (oriented) sum of the complex phases of the bond variables 
around a plaquette defines a U{1) flux. This flux is related to the solid angle formed 
by the spins and it vanishes in coUinear as well as in coplanar states. In such cases the 
phases can be therefore be gauged away and the (Atj) can be made real. For this reason 
complex bond variables are usually not observed . F^l 



February 2, 2008 4:30 



WSPC/Trim Size: 9in x 6in for Review Volume 



a 



Two-dimensional quantum antiferromagnets 27 

defined by bonds where the exchange Jij ^ is bipartite. This 
can also be true if the original lattice is a non-bipartite lattice 
but some bonds have Qij = so that the remaining lattice is 
bipartite. This is JJie case, for instance, in the J1-J2 model on 



the square lattice,lii-l in some regions of the phase diagram of the 
Shastry^utherland modcEj and in on the checkerboard Heisenberg 
model.E2fl Such configurations of the bond variables give collinear 
spin structures: spin-spin correlations can either be long-ranged 
(large k, Neel phase) or short-ranged but in both cases the magnetic 
structure factor is peaked at a simple wave- vector kg such that 2ko 
is a reciprocal lattice vector (feg — (tt, 0), ko ~ (0, tt) or fco — (t^, tt) 
in square geometries), 
ii) The lattice made of the sites connected by non-zero Qij bonds is 
not bipartite. This ha ppe ns in some phases of the J1-J2-J3 m odel 
on the square lattice,E£2l on the triangular or kagome latticesE^° 



in the Shastry-Sutherland model for some values of the exc hange 



parametersEj and on an anisotropic triangular lattice.c^ll Such 
mean-field states generically have planar but non-coUinear spin- 
spin correlations. 

It is simple to check that case i) preserve a global continuous symmetry 
while such a symmetry is absent in ii). Consider the following global gauge 
transformation in case i) : 

b^a 6..e'^ bj, -> bj^e-"^ (26) 

Qii' ^ Qii'^ Qjj' ^ Qjj'^ *^ (^'^) 

Qij Qij (28) 

where belongs to sublattice A and j,j' belongs to sublattice B. This 
transformation does not change the mean-field parameters Qij. The only 
fields affected by this transformation are those connecting two sites on the 
same sublattice.' They have a zero expectation values in case i) (or do not 
even exist if the physical lattice is itself bipartite) . For this reason it is pos- 
sible to make low-energy (and long-wavelength) gapless gauge excitations 
about the saddle-point by replacing the global staggered phase shift cj) of 
Eq. by a slowly varying (staggered) 0ij (Eq. |2^) . A gradient expansion 

'in a gauge theory language the fields Qui (resp. Qjj > ) of Eq. ^ transform like a charge-2 
scalar for the U{1) gauge field. Instead, from Eq. Eq, the bosons (spinons) carry a charge 
1. 
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of the effective action performed at the appropriate points in the BriUouin 
zone for the phase fluctuations only involves gradients of 9. The correspond- 
ing action is that of U{1) lattice gauge theory coupled to charge- 1 boson 
(spinons) . 

4.3.4. Topological effects - instantons and spontaneous dimerization 

So far only small fluctuations around the saddle point were considered and 
the contribution of topologically non-trivial gauge-field configurations were 
ignored. Consequently, the magnitude of the "spin" nc/2 was a continu- 
ous parameter and the information about the integer or half-integer (for 
instance) character of S as disappeared. Fmm Haldane's work on quantum 
spin chains and non-linear sigma modelsi^ it is well known that Berry 
phases in spin systems give rise to topological terms in the low-energy ef- 
fective action which can play a crucial role depending on the parity of 25". 

In non- linear sigma models in 2-1-1 dimensions the Berry phase van- 
ishes f or con figurations which are smooth on the scale of the lattice 
spacingE3E^ (unlike the 1-1-1 dimensional case). However "hedgehog" 
space-time singularitieal23 give non-trivial Berry phases. Read and Sachdev 
found that the closely related instantons of the effective U{1) gauge theory 
described before™ also play a crucial role in the physics of the Sp{N) (as 
well as SU{N)) spin models. 

The Berry phase associated to such a singular configuration depends on 
details of the lattice geometry. In the short-range ordered (tt, tt) phase of 
the square lattice antiferromagnet this Berry phase is a multiple of incTr/2. 
Although dealing with a gas of interactin g ( 1/r Coulomb-like potential) 
instantons is a difficult problem (see Ref. E£l and references therein), we 
can guess that the physics will depend on Uc mod 4. A detailed analysisllii 
shows that when Uc ^ mod 4 the instantons condense and spontaneously 
break the lattice translation symmetry. This generates a static electric field 
for the U{1) gauge field. Since the electric field is coupled to the difference 
of amplitudes of the bond variables, such state acquires spatially inhomoge- 
neous expectation values of the bond variables, it is a VBC and spinons are 
confined in pairs. In the J1-J2 model around J2/J1 — 0.5 the mean-field 
state is short-range ordered with correlations peaked at (tt, tt). A columnar 
dimerized state is predicted by Read and Sachdev from this analysis of the 



™In that gauge theory associated to the phases of the link variables an instanton corre- 
sponds to a tunneling event during which the total magnetic field piercing the lattice is 
changed by ±271. 
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fluctuations, in agreement with a number of numerical works on the SU{2) 
J1-J2 spin-i model. 



In a recent work by Harada, Kawashima and TroyerE2£l the phase dia- 
gram of the (first neighbor - unfrustrated) SU (N) antiferromagnet on the 
square lattice with ric = 1 was found to be in complete agreement with 
Read and Sachdev's predictions. They showed by quantum Monte Carlo 
simulations that for < 4 the systems is Neel ordered whereas it is a 
columnar VBC for N > 5. This provides an additional support to the field 
theory arguments described above. It also underlines that the mechanism 
of spontaneous symmetry breaking and formation of a VBC may come 
from quantum fluctuations only and that frustration is not always required 
(although it clearly enhances quantum fluctuations). 

On the other hand when ric = mod 4 the analysis of Read and Sachdev 
shows that fluctuations should not bring any broken lattice symmetry. 
Spinous are also confined and this state_closely resembles the valence-bond 
solid (VBS) proposed by Affleck et alM as a possible ground-state when 
the spin S matches the coordination number z according to 25 = mod z 
(see Q. 



4.3.5. Deconfined phases 

Now we suppose that, starting from a mean-field solution with collinear 
correlations (case i), a parameter of the original spin model is varied so 
that the mean field solution is changed and some bonds Qu' {i and i' be- 
long to the same sublattice) acquire a non zero expectation value (case ii) . 
In the framework of square lattice_antiferromagnets, a finite third-neighbor 
coupling (J3) would be needed.ll23 From the point of view of the long- 
wavelength gauge fluctuations (related to the continuum limit of the phases 
Oij) discussed above, the appearance of Qu' 7^ is equivalent to the con- 
densation of a (Higgs) boson of charge 2. This is a spontaneous break down 
of the global U (1) staggered symmetry of Eqs. p6|-p8|down to a Z2 one since 
the field Qu is not invariant iinr er Eq. p7| except if = or tt. Based on re- 
sults of Fradkin and Shenkert2£l concerning confinement in compact lattice 
gauge theories coupled to matter, Read and Sachdev argued that this Higgs 
mechanism suppresses the low-energy gauge fiuctuations and liberate the 
spinous. This confinement transition is described by a Z2 gauge theory. The 
suppression of the U{1) gauge fluctuations also forbids the condensation of 
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instantons discussed above and the ground-state remains unique" and bond 
variables have uniform expectation values. It is a genuine SL without any 
broken symmetry and deconfined spinons. 

5. Quantum Dimer Models 

In a previous section (§||) we showed that pairing spins- ^ into singlets at 
short distances is a rather natural way to overcome frustration in Heisen- 
berg antiferromagnets. QDM are defined in the Hilbert space of nearest- 
neighbor valence-bond (or dimer) coverings of the lattice. By construction 
these models focus on the dynamics in the singlet space and ignore mag- 
netic (gapped magnons or gapped spinons) excitations. For this reason they 
are (a priori) not appropriate to describe the physics of spin systems where 
magnetic excitations are gapless. 

The Hamiltonian of a QDM usually contains kinetic as well as potential 
energy terms for these dimers. Such Hamiltonians can often be simpler than 
their spin parents and are amenable to several analyti c treatrn ents because 



of their close relatio ns to c l a,ss ic3,l dimer problems ,[!£3t£al!±3 Ising models 



and Z2 ga.uge theory£lJ'tid'E13 T hese rn odels can offer simple descriptions 



of VBQilj as well as RVB liquids. tlSliiil It is in particular possible to write 



down some QDM that have a simple and exact VBC ground-state with 
spontaneous brok en symmetries (such as Rokhsar and Kivelson's model on 



141 



the square latticetd^ with attractive potential energy only - in which case 
the exact ground-state is very simple). Simple sol vabl e QDM which have a 
dimer-liquid ground-state can also be constructed, 

5.1. Hamiltonian 



The first QDM was introduced by Rokhsar and Kivelson.Elj It is defined by 
an Hamiltonian acting in the Hilbert space of first-neighbor dimer (valence- 
bonds) coverings of the square lattice and reads: 

^= E [-^(lii>(«l+H.c.)+y(in)(iii + ir:)(r:l)] (29) 

Plaqucttc 

Flipping two parallel dimers around a square plaquette is the simplest 
dimer move on the square lattice and the J terms precisely represent such 
dynamics. The V terms are diagonal in the dimer basis and account for an 
attraction or repulsion between nearest-neighbor dimers. These are the two 



"Except for a discrete topofogical degeneracy. 
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most local terms (respecting all lattice symmetries) which can be consid- 
ered.^ 

5.2. Relation with spin-^ models 

Sercnt interesting mappings between frustrated Ising models 



There exists. 

and QDMhBHi Here, however, we focus on the relations between QDM 
and SU{N) (or Sp{N)) spin models in which dimers are related to singlet 
valence-bonds. 

Overlap expansion. — A valence-bond state (product of two-spin sin- 
glets - belongs to the spin-^ Hilbert space) can be associated to any dimer 
covering.*^ Two such valence-bond stat es \a) and \b) are not orthogonal 
but, as first discussed by Sutherlandlilj their overlap decays exponentially 



with the length L of the loops of their transition graphs (defined in g5.3.l| ) 



as I {a\b) I = 2^"-^/^. Rokhsar and Kivelsonlilj introduced a formal expan- 
sion parameter x and replaced | (a|6) | by 2x^ . Their idea is that although 
X = for physical SU{2) spins, the physics of some models may be cap- 
tured by the first orders of a small x expansion. Truncating this expansion 
to order x" gives an effective Hamiltonian which contains local terms in- 
volving at most n dimers.*^ In this approach the dimer states of the QDM 
are in one-to-one correspondence with orthogonalized valence-bonds in the 
spin Hilbert space. 

Fluctuations about large-N saddle points. — From the argument above 
it could seem that the connection between spin-i models and QDM relies 
on a variational approximation: the spin Hilbert space is restricted to the 
nearest-neighbor valence-bond subspace. This connection is in fact prob- 
ably deeper, as some theories describing fluctuations about some large- iV 



^Thevwere originally derivediiJ as the lowest order terms of a formal overlap expansion 
(see § p.2| below) of the Heisenberg. In that calculation J ~ a:** and V . Notice that 
a three-dimer kinetic term (extending over two neighboring plaquettes) is generated at 
order and is not included in Eq. GSl. 

'^There is however a sign ambiguity (a valence-bond is antisymmetric under the exchange 
of both spins) that can be fixed by choosing an orientation on every bond. 
•^The signs of the non-diagonal (kinetic) terms of the effective QDM obtained by such 
small X expansion depends on the sign convention which was chosen to map valence- 
bonds to dimers. An important question is to know whether, at least at the lowest 
non-trivial order, a sign convention givine.4he same sign for all kin»ti« terms exists (as 
in Eq. ^). This is the case on the squar J-*^^ and triangular lattices 

"^This implicitly assumes that the valence-bond states are linearly independent. This can 
be demonstrated on the square lattice and appears to be the case on the triangular and 
kagome lattices for large enough sizes. The operator which orthogonalizes the valence- 
bond basis into the dimer basis is Q"^/^ where Q,^ i, = {a\b) is the overlap matrix. 
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saddle points corrections) are equivalent to (generalized) QDM. This 

mapping was discussed by Read and SachdevllHj for representation with 
m = 1 (number of rows in the Young tableau of the SU{N) representa- 
tion), N oo and 0{1), it leads to a generalized QDM where ric 
dimers emanate from each site. A QDM also describes 1/N corrections in 



the case of the fermionic'' SU{N) generalization 
model: 



of the Heisenberg 



= -^^Jy 4i?y+cst 



N 



where Bi. 



Ei 

<T=1 



(30) 
(31) 



and where the Cia are N flavors of fermions satisfying a constraint similar 
to Eq. |l| : 



N 

E 



N/2 



(32) 



Rokhsar showedE£fl that in the iV — > cx) limit "dimerized states" (or Peierls 
states) becomes exact ground-states of Eq. ^ for a large class of models.^ 
Quite naturally, 1 /N corrections will induce a dynamics into this subspace 
of dimerized states; it can be described by a QDM (with kinetic energy 
terms only at this order). At lowest order, on the square lattice, a kinetic 
term identical to the J term in Eq. ^ is generated and naturally favors a 
columnar or resonating-plaquette crystal (in agreement with a number of 
works on the spin-i model). The same arguments were discussed for the 
kagome lattice. tl3 In that case the leading corrections to the fermionic 



saddle point generate three-dimer moves around hexagons and stabilize a 



^ Young tableau with Uc = 1 column and m = N/2 rows. 

'Let Jo be the largest value of the Jij . Each dimerization (no site is left empty) where 
only bonds where Jij = Jo are occupied is a ground-state. Here a "dimer" between to 
neighbors i and j consists of a SU (N) singlet made with N fermions (one of each flavor) 

hoping back and forth between i and j. It is constructed from H^Li (j^la + '^jtr) l*^)- 
projecting out the components which do not satisfy Eq. Notice however that in the 
N ^ oo limit the relative fluctuations of the total number of fermion on each site are of 
order 1/V~N and can be neglected. 
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crystal of resonating.^ hexagons.^ This formalism was also applied to the 
checkerboard modelEa to conclude to the presence of a VBC phase. 



5.3. Square lattice 



The phase diagram of the Rokhsar and Kivelson's square lattice QDM is 
shown Fig. |^. Since a change in the signs of the basis dimer configurations 
can change J into — J (see Ref. 1113) we will choose J > without lost of 
generality. 





V/J 
— ► 



—0.2 ? 



Fig. 3. Schematic phase diagram of the the square lattice QDM. The popeifele location 
of an intermediate plaquette phase is taken from the work of Leung et al£^ 



5.3.1. Transition graphs and topological sectors 

In order to understand the particularities of QDM on bipartite lattices it 
is useful to describe how the set of dimer coverings splits into topological 
sectors. To do so we first have to define transition graphs: the transition 
graph of two dimer coverings c and c' is obtained by superimposing c and c' 
on the top of each other; it defines a set of non-intersecting loops covering 
the lattice. On each bond where the dimers of c and c' match, a trivial 
loop of length 2 is obtained. When the lattice is bipartite (two sublattices 
A and B) these loops can be oriented in the following way: any dimer 
belonging to c is oriented from A to B and dimers of c' are oriented B —^ 
A. The transition graph is then made of loops of the type A ^ B 
A B ■ ■ ■ . With periodic boundary conditions two winding numberall; 
—L/2 < rtxj^y < L/2 are associated to such a transition graph {L x L 
sites), fix (resp. Qy) is the net number of topologically non-trivial loops 



^To our knowledge however there is no clear evidence of such ordering in the spin-^ case. 
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(clockwise minus counterclockwise) encircling the torus in the x (resp. y) 
direction. 

Dimer coverings can be grouped into different topological sectors. By 
definition two dimer coverings belong to the same sector if they can be 
transformed into each other by repeated actions of local dimer moves (the 
transition graph associated to each movement does not wind around the 
whole system if it has periodic boundary conditions) . On the square lattice 
two-dimer moves are sufficient to connect any two states in the same sector; 
that is the Hamiltonian Eq. ^ is ergodic within each topological sector. 
In a torus geometry, c and c' belongs to the same sector if and only if 
their transition graph has winding numbers U.^ = = 0. The different 
topological sectors can be labeled by their winding numbers with respect 
to some reference columnar configuration. Their number is of order O(L^) 
for a system of linear size L. 

5.3.2. Staggered VBC for V/J > 1 

When V is sufficiently large the system tries to minimize the number of 
parallel dimers. The staggered configuration shown Fig. |^ has no such flip- 
pable plaquette. It is always a zero-energy eigenstate of Eq. ^ and becomes 
a ground-state for J > y. It breaks several lattice symmetries (four- fold 
degenerate) and is a VBC. 

The expectation value of the energy per plaquette satisfies min(0, V — 
J) < Eo/Np < max(0, V + J). For V/J > 1 this gives < Eo/Np and any 
zero-energy state saturates this lower bound and is a therefore a ground- 
state. One should however notice that it is possible to make zero-energy 
domain walls in this VBC since the state shown Fig. ^ is also annihilated by 
the Hamiltonian. No local dimer movement can take place in the staggered 
VBC (with or without domain walls). Each of these states form a topological 
sector with a single configuration which has jOj:! -|- = L/2 with respect 
to a columnar state. 

5.3.3. Columnar crystal for < 

When parallel dimers attract each other the system tries to maximize the 
number of flippable plaquettes. Columnar configurations as shown on the 
left side of Fig. ^ do maximize this number. Such a VBC is exactly realized 
for < and J — 0. Elementary excitations consist of a pair of (say) 
vertical dimers in a background of vertical columns of horizontal dimers. 
Such excitations are gapped {AE = 2\V\) and this VBC phase will survive 
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Fig. 4. Domain wall in a staggered VBC. 



to the inclusion of a finite J term. Notice that unlike the staggered VBC 
presented in the previous paragraph the columnar dimer configuration is 
not an exact eigenstate when J ^ Q. The exact diagonalizations data of 
Leung et allli suggest that the columnar phase may disappear around a 
critical value V/J ~ —0.2. 



5.3.4. Plaquette phase 

When the kinetic energy dominates {\V\ <^ J) the system will try to max- 
imize the number of resonating plaquettes |||) + |=). This can be achieved 
through the resonating plaquette crystal shown Fig. ^ and th e nu merical 
work (exact diagonalizations up to 8 x 8 sites) of Leung et alil3 suggest 
that this phase is realized in an interval —0.2 < V/J < 1. Although this 
model would not suffer from the fermion sign problem we are not aware of 
any QMC simulation on this model. 



5.3.5. Rokhsar-Kivelson point 

The point J — V (Rokhsar-Kivelson_(TlK) point) plays a special role. As 
remarked by Rokhsar and Kivelsonlll3 the equal-amplitude superposition 
of all dimerizations in a given topological sector is an exact ground-state. 
When J = V the Hamiltonian can be written as a sum of projectors: 

Hj=y=i =5]|vl',)(vl/p| (33) 
P 

\%)-\u)-\:z) (34) 

The linear superposition of all dimer coverings belonging to a given sector 

n 

|0) = ^ |c) (35) 
cen 
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is annihilated by Eq. ^ and is therefore a ground-state. The argument is 
the following. Consider a plaquette p and a configuration |c). If |c) has one 
or no dimer at all on the edges of p we have (5'p|c) = 0. If two dimers are 
present, then there exists a configuration |c') in the same sector which only 
differ from |c) by a two-dimer flip on p. In such a case the combination 
|c) + |c') is again orthogonal to |^'p). This shows that H |0) = 0. 

When open boundary conditions are considered (this restricts the topo- 
logical sector to Vl^ — Vty = 0) the RK state is the linear combination 
of an exponential number of configurations.'^ This is very different from 
the crystalline states considered so far where some periodic configurations 
were favored and it closely resembles Anderson's RVB picture. As we shall 
see this RK state is not a "true" liquid on the square lattice since dimer- 
dimer correlations are not short-ranged but algebraically decaying (~ 1/?'^) 
with distance. The calculation of dimer-dimer correlations in the R K state 



137 



139 
143 



(Eq- P.5D maps onto a classical dimer problem solved by Kasteleyn 
FisherE^a and Fisher and Stephenson.El^ From this Rokhsar and Kivelson 
constructed gapless excitations (in single-mode approximation) which dis- 
persion relation vanishes as at small mmaentum (the origin is taken 



at (tTjTt)). Quoting Rokhsar and Kivelsonjlia these excitations (dubbed 
"resonons" ) are the Goldstone mode of the gauge symmetry which allows 
the phases of the different topological sectors to be varied without changing 
the energy. Another mode of gapless excitations (around (tt, 0) and (0, tt) 
in the Brillouin zone), specific to the fact that the ground-state_has critical 
(algebraic) dimer-dimer correlations, was recently discovered.t^ 

The QDM on the square lattice is thus believed to be ordered (VBC) 
everywhere except at the RK point (J — V) where it has quasi long-ranged 
(critical) dimer-dimer correlations. 

5.4. Hexagonal lattice 

The QDM on the honeycomb lattice was studied by Moessner, Sondhi and 
Chandra,E£il in particular with Monte Carlo simulations. The phase diagram 
is very similar to the square lattice-case discussed above. It possesses three 
crystalline phases (Fig. |^) and has algebraically decaying dimer-dimer corre- 
lations at the Rokhsar Kivelson point (where the ground-state in each sector 



has the same form as Eq. 35). The absence of liquid phase (with exponen- 



tially decaying 2n-mer-2ri-mer correlations) is believed to prevail in bipar- 



''This is also true for periodic boundary conditions provided the two winding numbers 
do not scale like the linear size L. 
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tite lattices. This relation between the absence of a deconfined dimer liquid 
phase ' and the bipartite character of the lattice as been discussed by sex. 



eral authoral^ and is related to the possibility of a height representation 
of dimer coverings when the lattice is bipartite.-' 
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(36) 



Fouet et aZ0 studied the spin-i J1-J2-J3 model on the hexagonal lat- 
tice by exact diagonalizations and found evidences of a staggered VBC of 



the type predicted for V/J > 1 by Moessner et aL[!£il in the QDM. Other 
phases (Neel ordered phase and a possible short-range RVB SL) are also 
present in the spin-i modeLEZI 



5.5. Triangular lattice 

The most local dimer Hamiltonian on the triangular lattice contains kinetic 
and potential two-dimer terms on each rhombus; it was studied by Moessner 



'Bipartiteness seems to forbid deconfinement but not short-ranged dimer-dimer corre- 
lations. The 4-8 lattice (squares and octogons) is an example where dimer-dimer corre- 
lations are short-ranged. We thank R. Moessner for pointing this to us. On this lattice 
the equal-amplitude superposition of all coverings would be similar to an explicit VBC 
wave-function (thus confining). Such situations are only possible when the number of 
sites is even in the unit cell. 

J Consider a bipartite lattice with coordination number z. For each dimer covering we 
can associate integers (heights) on the dual lattice by the following rule. Set the height 
to be zero on a plaquette at the origin. The height is then defined on the whole lattice 
by turning clockwise (resp. counterclockwise) around sites of the A-sublattice (resp. B- 
sublattice) and changing the height by z — 1 when crossing a dimer and by —1 when cross- 
ing an empty bond. It is simple to check the difference of heights &h(x) = h\(x) — h2{x) 
between two dimerizations is constant inside each loop of their transition graph. Notice 
that the loops of a transition graph can be naturally oriented on a bipartite lattice. Then, 
&h{x) changes by +2 (resp. —2) when crossing a clockwise (resp. counterclockwise) loop 
of the transition graph. Columnar dimerizations have an averaged height which is flat 
and staggered ones have the maximum tilt. The winding numbers (Oa;,Oy) correspond 
to the averageheightdifference between both sides of the sample. The kinetic energy 
terms of Eqs. Ea and Bq change the height of the corresponding plaquette by ±2 and the 
potential terms {V > 0) favor tilted configurations. 
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Fig. 5. Phase diagram of the hexagonal QDM obtained by Moessner et aZii£il Although 
the VBC shown for V < do not have all dimers parallel to the same direction it is 
equivalent to the columnar VBC found in the square lattice QDM because it maximizes 
the number of flippable plaqpettes. It also corresponds to the ordering pattern predicted 
in the large- A'^ theory of Ref. As for the VBC obtained for V/ J > 1, it is the hexagonal 
counterpart of the staggered VBC of the square lattice (no flippable plaquette, exact 
eigenstate and maximum tilt in a height representation) . 



and Sondhi 



144 



r 



H.c. 



(37) 



where the sums run over all rhombi r of the lattice (with three possible 
orientations). This model was shown to possess (at least) three crystalline 
phases, including a columnar and staggered one (Fig. ^ as in the two pre- 
vious examples. An additional VBC (with resonating diamonds plaquettes) 
with a large unit cell (12 sites) was also predicted around V = 0. When 
V < and J = the ground-state is highly degenerate since it is possible, 
from an ordered columnar configuration, to shift all the dimers along any 
straight line without changing the number of flippable plaquettes (contrary 
to the square lattice case). However, an infinitesimal J is expected to lift 
this degeneracy and to order the ground-state in a columnar way. 

The phenomenology of t hese ordered phases is that of usual VBC and 
we refer to the original papeillU for details. 



5.5.1. RVB liquid at the RK point 

The new physics of this model appears through the existence of a liq- 
uid phase in the interval 0.7 V/J < 1. As for the two previous QDM 
the ground-states are exactly known at the RK point J = V. As before 
dimer-dimer correlations are obtained exactly at this point by a Pfaffian 
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calculation 



144 154 



but the result shows a fipMf correlation length. From 
their Monte Carlo simulations Moessner et a/.Edil argued that the spectrum 
is gapped at the RK point and that this gap persists down to V/J~ 2/3, 



that is over a finite range of coupling. This pic 
exact diagonalizations performed on this model. 



urn is consistent with the 



156 



157 



5.5.2. Topological order 

When the lattice is not bipartite the loops of a transition graph can no 
longer be oriented. The winding numbers ^Ix and are now defined as the 
(positive) number of non-trivial loop around x and y in the transition graph 
with a reference configuration (say columnar). These two integers are not 
conserved by local dimer moves. They are however conserved modulo two, 
which leaves only four sectors.'' Consequently the dimer liquid ground-state 
is four-fold degenerate at the RK point. This degeneracy holds exactly at 
the RK point even on finite-size samples but it is expected to hold in the 
thermodynamic limit in the whole liquid phase (0.7 ^V/J< 1). 

Conventional orders are often associated to a spontaneously broken sym- 
metry and lead to ground-state degeneracies in the thermodynamic limit. 
The four-fold degeneracy disc usse d above is the signature of some kind of 



order, called topological order .E3 The peculiarity of this order is that it is 
not associated to any local order parameter: a local observable cannot de- 
cide whether a given dimerization is in one sector or another. The existence 
of topological order is intimately associated to the fractionalized nature of 
the elementary excitations (see below). In the case of a RVB dimer liquid 

''The two-dimer moves included in Eq. ^ are not sufficient to guaranty ergodicity with 
each of the four sectors. Staggered states (12-fold degenerate - not 0{L) like on the 
square lattice) hajiiie-|no flippable plaquette but can be connected to other states with 
four-dimer moves 
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these excitations have been known to h e TRin.g vortices for a long timeE££[E£E 
(dubbed visons in the recent Hteraturell^ll^) . We will now discuss these 
excitations in more details in the framework of a QDM which realizes the 
same dimer liquid phase but for which not only the ground-state but all 
the eigenstates are known exactly. 



5.6. Solvable QDM on the kagome lattice 



An exactly solvable QDM on the kagome lattice was introduced by D. Ser- 



ban, V. Pasquier and one of usiifil. It offers a very simple and explicit 



realization of the ideas discussed above (visons, topological order etc. 



5.6.1. Hamiltc 



The kagome lattice QDM introduced in Ref. IHil contains only kinetic terms 
and has no external parameter. The Hamiltonian reads: 

n^-Y,a^{h) (38) 



32 



where ct^(/i) = Idc^W) {dc.ih)\ + H.c (39) 

a=l 

The sum runs over the 32 loops on the lattice which enclose a single hexagon 
and around which dimers can be moved (see Tabic |l] for the 8 inequivalent 
loops). The shortest loop is the hexagon itself, it involves 3 dimers. 4, 5 
and 6-dimers moves are also possible by including 2, 4 and 6 additional 
triangles (the loop length must be even). The largest loop is the star. For 
each loop a we associate the two ways dimers can be placed along that 
loop: \da{h)) and \da{h)y Notice that a^{h) measures the relative phases 
of dimer configurations displaying respectively the da{h) and da{h) patterns 
in the wave function. 



5.6.2. RK ground-state 

As for the QDM discussed previously the ground-state of this Hamiltonian 
is the equal amplitude superposition of all dimer coverings belonging to a 
given topological sector (as on the triangular lattice there are four sectors). 
This can be readily shown by writing 7i as a sum of projectors: 

32 
h a=l 
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Table 1. The 8 different classes of loops which can 
surround an hexagon of the kagome lattice. Including all 
possible symmetries we find 32 possible loops. The first 
column indicates the number of dimers involved in the 
coherent motion around the hexagon. 
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, >::< 
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where is the number of hexagons on the lattice. When expanding the 
products the diagonal terms give a simple constant since 



32 

E 

a=l 



da) {da\ + \da) {da\ 



1 



(41) 



This reflects the fact that, for any dimerization, the dimers on hexagon h 
match one and only one of the 2 x 32 patterns {da, da} ■ 

Unlike the square or triangular case, the RK ground-states |0) = 
Ecsn I'') ^^'^ degenerate with some staggered VBC' This means that 
the Hamiltonian of Eq. ^is not at a phase transition to a VBC. As we will 
explain it is inside a liquid RVB phase. 

The RK wave-function can be viewed as dimer condensate. It is similar 
to the ground-state of liquid ^He which has the same positive amplitude 
for any configuration and its permuted images.llS3 An important difference, 
however, is that the QDM state is incompressible and cannot sustain acous- 
tic phonons. This can be related to the fact that the U{1) symmetry of the 
Bose liquid is absent in the QDM on non-bipartite lattices. It is replaced 



'Because resonances loops of Itmgj^h up to 12 are present the dynamics is ergodic in each 
of the four topological sectors. |Hj 
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instead by a discrete Z2 gauge symmetry (see §5.6.7 below) 



5.6.3. Ising pseudo-spin variables 

The kinetic energy operators defined in Eq. ^ commute with each other. 
This is obvious when two such operators act on remote hexagons but it also 
holds for neighboring ones. This property can easily be demonstrated with 
the help o f th e arrow representation of dimer coverings introduced by Zeng 
and Elser.E£3 This mapping of kagome dimerizations to arrows on the bonds 
of the honeycomb lattice is illustrated Fig. ^ Each arrow has two possible 
directions: it points toward the interior of one of the two neighboring tri- 
angles. If site i belongs to a dimer {i,j) its arrow must point toward the 
triangle the site j belongs to. A dimer covering can be constructed from 
any arrow configuration provided that the number of outgoing arrows is 
odd (1 or 3) on every triangle. 



Fig. 7. A dimer covering on the kagome lattice and the corresponding arrows. Dashed 
lines: honeycomb lattice. 



The operators have a particularly simple meaning in terms of the 
arrow degrees of freedom: a'^lh) flips the 6 arrows sitting on /i.™ It is then 
clear that the commute and that a^{Kf' = 1. In fact these operators can 
be used as Ising pseudo-spin variables and the Hamiltonian now describes 
non-interacting pseudo-spins in a uniform magnetic field pointing in the x 
direction. In the ground-state we have a^{K) = 1 on every hexagon. 

■"Flipping all the arrows around any closed loop (such as around an hexagon) preserves 
the local constraint imposed on arrow configurations. Flipping the arrows around a 
topologically non-trivial loop changes the topological sector. 
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5.6.4. Dimer-dimer correlations 

The ground-state is the most possible disordered dimer hquid as the dimer- 
dimer correlations strictly vanish beyond a few lattice spacings. Such cor- 
relations can be computed by the Pfaffian method. On the kagome lat- 
tice the determinant of the Kasteleyn matrix (which is directly related to 
the partition fu nctio n of the classical dimers problem) is exactly constant 



in Fourier space.ll£j Since dimer-dimer correlations are obtained from the 
Fourier transform of the inverse of this determinant, they turn out to be 
strictly zero beyond a few l atti ce spacings (as soon as the two bonds do not 



touch a comr non tri angle). E-jj This result can also be obtained by a sim- 
pler argumentlilil'E£3 using the cr^ operators. This result is related to the 
kagome geometry.'^ This absence of long-ranged dimer-dimer correlations 
demonstrates that the RK state is a dimer liquid and that it breaks no 
lattice symmetry. 

On the triangular lattice, even at high temperature, dimer-dimer cor- 
relations decay exponentially with distance but these correlations remain 
finite at any distance. On the square lattice such correlations are even larger 
because they decay only as a power law. This means that the infinite hard- 
core dimer repulsion makes QDM non-trivial even at infinite temperature; 
dimers cannot be free when they are fully-packed. From this point of view 
we see that the kagome lattice is particular: it is as close as possible to a 
free dimer gas, except for non-trivial correlations over a few lattice spac- 
ings. This is a reason why dimer coverings on the kagome lattice can be 
handled with independent pseudo-spin variables and why the RK state on 
this lattice is the most possible disordered RVB liquid. 



5.6.5. Visons excitations 

The cr^ operators can be simultaneously diagonalized but they must satisfy 
the global constraint Y\hCF''^{h) = 1 since this product flips every arrow 
twice. It must therefore leave all dimerizations unchanged. The lowest ex- 
citations have therefore an energy 4 above the ground-state and they are 
made of a pair of hexagons a and 6 in a (T^(a) =jxl(i) = ^1 state, a and h 



are the locations of two Ising vortices (or wisoni£lrE£3) . As remarked before 
this means that the relative phases of the configurations with da{h) and 
da{h) patterns have now changed sign. The corresponding wave- function is 

"The m(^de^ of Eg. can be generalized to any lattice made of corner-sharing 
triangles. tlil 
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obtained in the following way. Consider a string which goes from a to 6 
(see Fig. ||) and let il(a, b) be the operator which measures the parity ±1 of 
the number of dimers crossing that string. ri(a, b) commutes with all (T^(ft,), 
except for the ends of the string: a^{a)^{a,h) — ~Q{a,b)a^ {a). A dimer 
move changes the sign of ri(a, 6) if and only if the associated loop crosses 
the string an odd number of times, which can only be done by surrounding 
one end of the string. This shows that SI (a, 6) flips the in a and b.° 
As the RK ground state |0), ri(a,5)|0) is a linear combination of all dimer 
configurations belonging to one sector. However the amplitudes are now 1 
and —1 depending on the number of dimers crossing Q. This wave- function 
therefore has nodes, it is an excited state of energy 4 with two vortices in a 
and b. It is easy to see that a different choice V.' for the string connecting a 
and b gives the same state up to a global sign which depends on the parity 
of the number of kagome sites enclosed by SI U SI'. 

These vortex excitations carry a Z2 charge since attempting to put two 
vortices on the same hexagon does not change the state. Such excitations 
are not local in terms of the dimer degrees of freedom. Indeed, determining 
the sign of a given dimerization in a state with two visons which are far 
apart requires the knowledge of the dimer locations along the whole string 
connecting the two vortex cores. In this model the visons appear to be static 
and non-interacting. This is a particularity of this solvable model but the 
existence of gapped vison excitations is believed to be a robust property of 
RVB liquids. In more realistic models the visons will acquire a dynamics 
and a dispersion relation but will remain gapped. p They will also have some 
interactions with each other but should remain deconfined. This property 
is particularly clear in the kagome QDM: visons are necessarily created by 
pairs but the energy is independent of their relative distances. 

The Ising vortices also offer a simple picture of the topological degen- 
eracy. Consider a ground-state |-|-) of the model which lives in the sector 
where the winding number f2j, (with respect to some arbitrary but fixed 
dimerization) is even. Another ground-state |— ) is obtained in the odd-f2y 
sector. Now consider the combination |0) = |+) + |— ) and apply the opera- 
tor fl(0, Lx) corresponding to a closed loop surrounding the torus in the x 
direction. This amounts to creating a pair of nearby visons at the origin, tak- 
ing one of them around the torus in the x direction and annihilating them. 



°Up to a global sign (reference dependent) Q{a, b) is equal to a" (a)a-^ (b) where the a" 
operators are those introduced by Zeng and Elser. 

Pit is possible to add potential energy terms to Eq. to drive the system outside of the 
liquid phase and this transition corresponds to a vison condensation. 
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Fig. 8. A pair of visons (located in a and b) is created by applying to the RK wave- 
function a factor (—1) for each dimer crossing the string Q. The dimerization shown 
there on the left appears in the linear superposition of the two-vison state with the sign 
— 1 whereas the one on the right has the sign +1. 



This can also be viewed as the creation of a vison in one hole of the torus 
(with no energy cost). It is simple to check that n{0, Lx)\0) = |+) — |— ) (up 
to an irrelevant global sign). This provides a simple relation between the 
vison-pair creation operator and the existence of two topologically distinct 
ground-states |-f) -I- |— ) and |+) — |— ). 



5.6.6. Spinous deconfinement 

We assume that dimers represent "dressed" singlet valence-bonds, as in the 
overlap expansion (§5.2). Since the Hilbert space is made of fully-packed 
dimer coverings the model of Eq. ^ only describes spin-singlet states. How- 
ever, as any QDM, it can be extended to include static holes or spinous. 
Configurations with unpaired sites (spinon or holon) are now allowed but 
the kinetic terms of the original Hamiltonian which loop passes on an empty 
site gives zem. Consider a system w ith t wo static spinons in x and y. As on 
the squarelil3 and triangular latticegll^ at the RK point the exact ground- 
state \x,y) remains the sum of all dimer coverings and the ground-state 
energy is independent of the distance between the two spinons (except at 
very short distance if they belong to a common hexagon) . This is a first in- 
dication that RVB spin liquid has deconfined spin-i excitations (spinons). 
In the QDM language these excitations are simply unpaired sites in a dimer 
liquid background. Such unpaired sites are necessarily created by pairs but 
they can then propagate freely (no attractive potential) when they are suf- 
ficiently far apart. 

Another calculation allows to test the deconfinement properties of a 
dimer liquid. We consider the state Itfj) — X^r^ol^''") where |0,r) is the 
(un-normalized) state with two spinons in and r. The probability to find a 
spinon in r in the lip) can be obtained by the relatively involved calculation 
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of the monomer correlation'^ withPfafRans. One the square lattice this 



probability goes to zero as l/y^.tla This shows that the second spinon 



is (quasi-) confined in the vicinity of the first one on the square lattice 
because escaping far away represents a large "entropy" cost in the dimei 
background. On the triangular lattice it goes exponentially to a constant-li^ 
This result is a signature of deconfinement. In fact the same signature can 
be obtained on the kagome lattice without any technic al c alculation since 



167 



the monomer correlation is exactly 1/4 at any distance. 

If unpaired sites are allowed one can describe spinous or holons. Unfor- 
tunately in the presence of simple kinetic energy terms for these objects the 
model can no longer be solved. However one can consider a static spinon and 
its interaction with visons: when the spinon is adiabatically taken around a 
vison the dimers are shifted along a path encircling the vison. Because the 
vison wave-function is particularly simple in this model it is easy to check 
that this multiplies the wave- function by a fac tor —1. This is the signature 



of a long-ranged statistical interactionE£3E£3 between visons and spinous 
(or holons). In more realistic models, as long as the visons are gapped ex- 
citations the spinous are expected to be deconfined. On the other hand if 
the visons condense their long-ranged statistical interaction with spinous 
frustrates their motion. This is no longer true if they propagate irupairs, 



in which case they are not sensitive any more to visons (see Ref. lilil for 
an extension of the present QDM with a vison condensation). This sim- 
ple physical picture illustrates the relation between vison condensation and 
spinon confinement. 



5.6.7. Z2 gauge theory 

The forces responsible for confinement are usually associated to gauge fields 
and their fluctuations. Where as r/(1 ) compact gauge theories are generically 
confining in 2 -I- 1 dimensions Z2 gauge theories are known to possess 

deconfined phases.li£3 For this reason some attention has been paid to the 
connexicma, between Z2 theories and fractionalized phases in 2D electronic 
systems. 



161 



It is knowi£13 that QDM can be obtained as special limits of Z2 gauge 
theories, the gauge variable being the dimer number on a bond. However, 
on the kagome lattice this connexion can be made exact and completely 



"JRatio of number of dimer coverings with two holes in and r to the number without 
hole. 
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myer 



14] 



In 



explicit since there is a one to one correspondence between dimer 
ings and physical states (i.e. gauge-invariant) of a Z2 gauge theory, 
this mapping the gauge fields are Ising variables living on the link of the 
honeycomb lattice {i.e. kagome sites) and are constructed from the ar- 
rows described previously. As for the constraints of gauge invariance they 
correspond to the odd parity of the number of outgoing arrows on every 
triangles. The tr^ operator used to define a solvable QDM translate into a 
gauge-invariant plaquette operator for the gauge degrees of freedom (prod- 
uct of the Ising gauge variables around an hexagon). With this mapping the 
visons appear to be vortices in the Z2 gauge field and the solvable model 
of Eq. |3^ maps to the deconfincd phase of the gauge theory. 

5.7. A QDM with an extensive ground-state entropy 

So far we have discussed QDM that realize either spontaneous VBC, critical 
states or RVB liquids. We wish to mention here that these three scenar- 
ios may not be the only possible ground-states for QDM. In particular, a 
QDM on the kagome lattice with an extensive ground-state entropy was re- 
cently discussed.ElJ This model was introduced from the observation that 
the dimer kinetic energy terms arising from an overlap expansion ( ^5.2| ) gen- 
erally have non trivial signs as soon as resonance loops of different lengths 
are considered. It was then realized that such signs (which make the QDM 
no longer appropriate for QMC simulations) can lead to qualitatively new 
phases diff erent from VBC Or RVB liquids. The Hamiltonian introduced 
in Ref. IlliJ is similar to that of Eqs. |38[]39| except that each resonance loop 
a is now included with a sign (—1)"° where = 3, • • • , 6 is the number 
of dimers involved: 

n = [\dc{h)) {da{h)\ + \d^{h)) {d^{h)\] (42) 

h 

These signs are precisely those appearing in the overlap expansion (at the 
order of one hexagon) of the Heisenberg model on the kagome lattice. This 



expansion was carried out by Zeng and Elserlllll in an insightful paper which 
laid the basis of the analysis of the kagome antiferromagnet in the first 
neighbor valence-bonds subspace. 

Although not exactly solvable, the Hamiltonian of Eq. ^ was shown 
to be a dimer liquid (short-ranged dimer-dimer correlation) and to have 
a huge ground-state degeneracy ~ 

2^/6 ^ 1.122^ {N is the number of 
kagome sites) . In additi on, severs,! numerical indications pointed to a critical 



behavior of this systemJiiiJ with a possible algebraic decay of energy-energy 
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correlations/ It was argued that the effective QDM describing the singlet 
dynamics of the spin-^ Heisenberg antiferromagnet on the kagome lattice 
could be close (in parameter space) to Eq. If correct, this sheds light 
on the very large density of singlet states observed at low energy in the 
numerical spectra of that spin model (see Sec. 0). 

6. Multiple-spin exchange models 

6.1. Physical realizations of multiple-spin interactions 
6.1.1. Nuclear magnetism of solid ^ He 

Solid ^He was the first mag netic syste m in which the importance of MSE 



interactions was recognized.blj'l H '"i Due to the large zero-point motion 



of the atoms about their mean positions, tunneling events during which 2, 
3 or 4 atoms exchange their positions in a cyclic way are frequent. These 
processes generate an effective interaction between the (nuclear) spins which 
can be written 

^ = E -JP^-^f {Psvrn + P^pl) (43) 
P 

where the sum runs over permutations P, Jp > is the exchange frequency 
of the associated tunneling process (in real space) and Pspin acts on the 
Hilbert space of spin-i located on the site of the crystal. The sign —(—1)^ 
depends of the signature of the permutation P and is a consequence of 
the Pauli principle. For a cyclic permutation involving n spins this sign is 
just (—1)" and is responsible for the ferromagnetic character of processes 
involving an odd number of spins. For spin-^ particles, two and three-spin 
exchange terms reduce to the familiar Heisenberg interaction: 

Pi2 = 25i • 5^2 + ^ (44) 

^123 + P32I = P12 + ^23 + P3I ^ 1 (45) 

but this is no longer true for n > 4: 

^'1234 + ^4321 = ^12^34 + -Pl4^23 ^ P13P24 + Pl3 + ^24 — 1 (46) 

which can be expressed (thanks to Eq. ^4[ ) as a sum of terms with two and 
four Pauli matrices. 



■^Notice that a one-dimensional analog of Eq. |42| can be defined and exactly maps onto 
the Ising chain in transverse field at its critical point. 
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■^He can form solid atomic mono-layers with a triangular geometry when 
adsorbed on a graphite substrate at ultra low temperatures (milli Kelyir 
range). This 2D magnet has been studied for a long time (see Refs. 



175 
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and references therein) and the importan ce nf M SF, interactions involving 
up to six atoms has now been recognized. tZZIIlZS The exchange frequencies 
of the most importa nt processes ^ have been computed by Path Integral 
Monte Carlo (P 
been carried out 



MO 



177 



79ll8qJl8Uh82 
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(analytic WKB calculations have also 
as a function of the density. The proposed MSE 
Hamiltonian describing the magnetic properties of this 2D quantum crystal 
reads 



n = {J2-2J,)Y,Pi2 + J4 (A...4 + H.C) (47) 

-J5 (A...5+H.c) + Je Y (^i...6 + H.c) 



where Eq. ^ was used to absorb the three-spin terms into an effective first- 
neighbor Heisenberg exchange = J2 — 2J3. At high density the hard- 
core potential between Helium atoms only leaves three-body exchanges 
possible {J3_^ Jn^a) and Eq. ^ reduces to a first neighbor Heisenberg 
ferromagnetlL^ as observed experimentally for the first time by Franco 
et aZ.E£3 in high-density layers. On the other hand the second layer so- 
lidifies at lower density anH higher order exchange terms cannot be ig- 
nored.*^ jEIMC simulationslH3 and high-temperature fits of the experimen- 
tal dataE-^ showed that the relative strength of two- and four-spin terms 
if roughly /J4 ^ —2 in the low-density second layer solid. The J2-J4 
model was studied by exact diagonalizations in this region of parameter 
space and evidences f or a sVi nrt- ranged RVB SL phase with no broken sym- 



metry were obtained. ll23'Ei3't£2l The ultra-low temperature measurements 
of specific heattSS and uniform susceptibilityfi^ are not incompatible with 
such a spin liquid phase but the spin gap, if any, has not yet been observed 
and should be smaller than 100/iK. 



"The first layer is then so dense that exchange is strongly suppressed. The first layer can 
also be replaced by an *He or HD mono-layer. 
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6.1.2. Wigner crystal 

The Wigner crystal is another fermionic solid with a triangular geometry 
where MSE interactions can play an important role. At very low density 
the Coulomb energy dominates, the crystal is almost classical and MSE 
interactions are very small. Exchange frequencies ./p can he co mputed in 
this regime by a semi-classical (WKB) approximatior£III^t22lll21l and, as for 
the high density solid of "^He, three-body exchanges dominate and give rise 
to ferromagnetism. However, at higher density and close to melting, PIMC 
calculations of the exchange frequencyll23 showed that the magnetism may 
be described by a MSE model with parameters ( J^^ and J4) close to those 
where the triangular MSE model is expected to be a RVB SL. Unlike the '^He 
case, the particles (electrons) are charged and an external magnetic field 
has also an orbital effect, it introduces complex phases in the exchange 
energies: P + P^^ — > e*"P -I- e~'^°'P~^ where the angle a = 27r0/0o is 
proportional to the magnetic flux (j) passing through the area enclosed by 
the exchange trajectory and S n is th e unit flux quantum. This can give 



with the Zeeman effect (see Ref. 
to explore this physics). 



rise to very rich phase diagramsE£jll£d where complex MSE terms compete 
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for some early experimental attempts 



6.1.3. Cuprates 

The possibility of signiflcant four-spin exchange around square Cu pla- 
quettes Jil copper oxide compounds was first suggested by Roger and 
Delrieu.t^il They interpreted the anomalously large width of Raman scat- 
tering spectra as a signature of four-spin exchange in this copper ox- 
ide superconductor. The importance of these MSE interactions in Cu02 
planes (J4 ^ O.25J2) has then been emphasized by a number of groups 
and in difF erent rnate rials and by different experimental and theoretical 
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oc 
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' ring exchange also plays a sig- 
'E23 For instance, exchange pa- 



approaches 1 

nificant role in ladder compounds 
rameters with values Jmng — Jieg = 110 meV and Jring = 16.5 meV were 
proposed f or La.RC aaCiiT/iOii based on the dispersion relation of magnetic 
excit at ions .E£3c23 



6.2. Two-leg ladders 

Numerous works were devoted to ladder models with four-spin interactions. 
These include general bi-quadratic interactions as well as models with ring- 
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exchange terms. We will only discuss here the simplest of these MSE models: 

n 

+ KY, (^234 + H.c) (48) 

□ 



Thanks to several studie^23E23'E2j'E£3 the phase diagram of this Hamilto- 



nian is now rather well understood and five different phases were identified. 

• Ferromagnetic phase. The ground-state is fully polarized. This 
phase includes the (J = —l,K = 0) and the {J = 0,K = — 1) 
points. 

• Rung-singlet phase. This phase includes the ground-state of the 
ladder without MSE term (J = 1, if = 0). The spectrum is gapped 
and the ground-state is iininue . A moderate K/J > 0.23 ± 0.03 



is 11 ni 

206 p07 



destroys this phasetfHJ^ElacHIIcHa in favor of the VBC below 



Staggered VBC with dimers on the legs. In one of the two de- 
generate ground-states the dimerized bonds are (2n, 1) — {2n+ 1 ,1) 
and {2n+l,2)- {2n + 2, 2). The VBC disappears for K/J > 0.5 



205 



Such a staggered VBC was first predicted in the framework ofa lad- 
der with bi-quadratic interaction by Nersesyan and_IsvelikO Us- 
ing Matrix-Product Ansatz, Kolezhuk and MikeskaEJ constructed 
models which are generalizations of Eq. ^ and which have exact 
ground-state with long-ranged staggered dimer correlations. In this 
phase the magnetic excitations are very different from the magnon 
excitations of the rung-singlet phase above. Here the excitations 
do not form well-defined quasi-particles but a continuum made of 
pairs of domain walls connecting two dimerized ground-states.@'0 
Scalar chirality phase. The order parameter is {Sn.i ■ {Sn,2 x 
'5'n+i.2)) and it spontaneously breaks the time-reversal symmetry 
and translation invariance. The ground-state is t wo- fold degenerate 



up to the next transi tion at K/J ~ 2.8 ± 0.3.Lr3 There exists a 



duality transformatioiiEHJ'E2j which maps the scalar chirality order 
parameter onto the dimer order parameter of the VBC above.' 
Applying such a transformation to the exact VBC ground-states 
mentioned above, models with an exactl y know n ground-state and 



scalar chirality LRO can be constructedE2J'c2a Although chiral SL 



have been much discussed in the literature, this is to our knowledge 

*The Hamiltonian of Eq. |48| is self-dual at 2K = J. 



February 2, 2008 4:30 



WSPC/Trim Size: 9in x 6in for Review Volume 



a 



52 



G. Misguich and C. Lhuillier 



the first realization of such a phase in a SU(2) symmetric spin-^ 
model. 

Short-ranged ordered phase with vector-chirality correla- 
tions. The strongest correlations are ((5'„^i x 8,1,2) ■ {Sn'.i x Sn'.2)) 
but they remain short-ranged. The spectrum is gapped and the 
ground-state is unique. This phase includes the pure K ~ 1 model 
where J = 0. This phase is related by the_diiality transformation 
discussed above to the rung-singlet phase.c2a This transformation 
indeed relates the Neel correlations {{Sn,i — Sn,2) • {Sn's — Sn\2)) 
(which are the strongest ones in the rung-singlet phase) to the 
vector-chirality correlations. Close to the ferromagnetic phase (J < 
0) one observes a crossover to a region where the strongest corre- 
lations are ferromagnetic spin-spi n co rrelations along the legs and 
antiferromagnetic along the rungs. 



6.3. MSE model on the square lattice 

The phase diagram of the Hamiltonian ^ o n th e two-dimensional square- 
lattice has been recently studied by LauchliE23 by exact diagonalizations. 
Neel, ferromagnetic, columnar VBC and staggered VBC phases were iden- 
tified, as in the ladder model above. In addition, a nematic phase character- 
ized by long-ranged vector chirality correlations (alternating spin currents) 
was found around the K = 1, J = point. To our knowledge this is the 
first microscopic realization of a nematic order in a two-dimensional spin-i 
model. 



6.4. RVB phase of the triangular J2 J4 MSE 

Because of its relevance to solid "^He films and Wigner crystals, 
MSE 



rondel 


on 


th 


^ trianmi 


sS^" 
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213 
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the 

. ar lattice has been the subject of many 
We will discuss here some properties of the 
simplest MSE model with up to four-spin cyclic exchange interactions 
(J2 — 2J3 and J4 only in Eq. 48). The classical phase di ag ram (Fig. ^ 
of this model has been studied by Kubo and c£>Ilaboratora2i3^Ei2l and the 



quantum one has been roughly scanned in Ref.li£3: we will mainly focus on 
the short-ranged RVB spin liquid (see Fig. ^), which is probably the first 
RVB SL encountered in an SU(2)-symmetric spin model. 
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Fig. 9. Classical (top) and quantum (bottom) phasesi-elnthe J2 — J4 MSB Hamiltonian. 
The classical model was studied by Kubo and MomoiElj and is based on a variational 
approach. The quantum phase dia groim m ti ne simplest scenario compatible with the 
exact diagonalizations data of Rofs. |i£3llij'E£j ■ While not completely understood, in the 
type-II spin-liquid region the spectra are characterized by a very large number of singlet 
excitations below the first triplet state. This is not the case in the RVB phase. 



6.4.1. Non-planar classical ground-states 

It is well-known that an Heisenberg model (with possible second-neighbors, 
third-neighbors, ... interactions) on a Bravais lattice always admits a planar 
helical ground-state at the classical level. This is no longer true when MSE 
are present and finding the classical ground-state for arbitrary J2 and J4 is 
an iiTiso]ve r| pr oblem A mean-field phase diagram was obtained for the clas- 
sical modelEij but very few exact results are known. In the neighborhood 
of J4 = 1 J2 = the classical ground-state is known to be a four-sublattice 
configurations with magnetizations pointing toward the vertices of a regu- 
lar tetrahedron.cii] This is a quite interesting model where the ground-state 
spontaneously breaks a discrete Ising symmetry associated to the sign of 
the triple product Si ■ {S2 x S3) around a triangle. This broken symmetry 
gives rise to a finite-temperaiiire phase transition which has been observed 



in Monte Carlo simulations. ciil This phenomena is similar to the transition 



predicted in the (tt, 0) phase of the J1-J2 model on the square latticed 



6.4.2. Absence of Neel LRO 

The classical ground-states at J4 = 1, ,/2 = are tetrahedral configurations. 
Although this phase a pp ears to be stable within the framework of linear 



exact 



spin-wave calculationaii-l or Schwinger-Boson mean-field theory 



diagonalizations indicate that the magnetic LRO is washed out by quantum 
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fluctuations-EHI'tifl Thexhiral order predicted to survive at long distances 
and finite temperaturesEill in the classical system for J2 = is also likely 



114 



to be washed out by quantum fluctuations. 

When J2 — I a relatively small amount of J4 0.1 is sufficient to de- 
stroy the thi£e-sublattice Neel LRO realized by the first-neighbor Heisen- 



berg modeLliSll The nature of the phase on the other side of this transition is 



not settled but the finite-size spectra display a large density of singlet ex£i 



187 



tations at low energy which could be reminiscent of the kagome situation 

From exact diagonalizations^imlp 36 sites) no sign of Neel LRO could 
be found at J4 = 1, J2 — 2.ll23lli3 In addition, the finite-size analysis 
showed that the spin-spin correlation length is quite short at J2 = — 2 and 
J4 = 1 and a spin gap of the order of A ~ 0.8 exists at this point. Much of 
the numerical effort to elucidate the nature of the MSE ground-state was 
concentrated on this point because it is close to the parameters realized in 
low-density '^He films (when higher order exchanges are neglected). 



6.4.3. Local singlet- singlet correlations - absence of lattice symmetry 
breaking 

Having excluded the possibility of a Neel ordered ground-state at J4 — 
1, J2 = —2 it is natural to look for a possible VBC. Because of the complex- 
ity of the MSE Hamiltonian it is not clear what kind of spatial order should 
be favored. From the analysis of dimer-dimer correlations (see Fig. |l^) it ap- 
pears that parallel valence-bonds repel each-other at short distance. This is 
similar to what is observed in the staggered phase of the J2-J4 MSE ladder 
and square-lattice models. For this reason it appears that a plausible VBC 
would be the staggered VBC encountered in the triangular QDM ior V > J 
(§ ^.5[ ). However this scenari o se ems difficult to reconcile with the weakness 
of dimer-dimer correlationsJiiJ In addition, the low-energy singlet states 
and their quantum numberscU do not reflect the 12-fold quasi-degeneracy 
that should be present if the system was to spontaneously break some lat- 
tice symmetry according to a staggered VBC pattern. Small systems usually 
favor ordered phases because low-energy and long-wavelength fluctuations 
that could destabilize an ordered state are reduced compared to larger sys- 
tems. From the fact that the finite-size spectra do not show the signatures 
of a staggered VBC symmetry breaking it is unlikely that the MSE model 
could develop a VBC of this kind in the thermodynamic limit. 
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Fig. 10. Dimer-dimer correlations in the ground-state oi^4lie J2-J4 MSE model on the 
triangular lattice (28 sites) at J2 = — 2, J4 = 1 (from Ref. Ilij). Numbers are proportional 
to (dodx) — {do){dj:) where the operator dj, projects onto the singlet space of the bond 
X and do refers to the reference bond (1,28). These results shows a clear tendency for 
repulsion between parallel dimers. 

6.4.4. Topological degeneracy and Lieb-SchuUz-Mattis Theorem 

Because no VBC phase could be identified in the MSE model at J4 = 
1, J2 = —2 the numerical data were compared with the predictions of an 
RVB liquid scenario. 



215 



In one dimension a famous theorem due to Lieb, Schultz and Mattis 
(LSM) states that in a one-dimensional spin system with an half-integer 
spin in the unit cell there is at least on excited state collapsing to the 
ground-state in the thermodynamic limit (periodic boundary conditions). 
There are in fact several arguments suggesting that this theorem migh t. 



at least to some extent, also apply to higher dimensions Ei3'tiZrll£ZJEL3'l£i3 If 
that is the case a gapped system with an odd integer spin in the unit cell 
must have a degenerate ground-state. The simplest scenario to explain this 
degeneracy is a translation symmetry breaking. One could think that this 
would rule out the possibility of any (translation invariant) RVB liquid in 
such models. This is incorrect because a ground-state degeneracy can have 
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a topological origin on a system with periodic boundary conditions, as we 
discussed in the framework of QDM (§|^). Such a phase is characterized 
by a four-fold topological ground-state degeneracy when the system is on 
a torus. That degeneracy allows the system to fulfill LSJVL's requirement 
without any spontaneous translation symmetry breaking. 
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On a finite-size system the topological degeneracy is only approximate 
but some constraints exist for the gua.n tum numbers (momentum in partic- 



ular) of the quasi-degenerate multiplet.E£ZI A system with periodic bound- 
ary conditions with an even number of sites but an odd number of rows is 
expected to have two ground-states with differ by a momentum tt in the 
direction parallel to the rows, in close analogy to the the LSM theorem in 
dimension one. The numerical spectra of the MSE model exhibit a set of 
three singlet energ y lev els collapsing onto the ground-state when the sys- 
tem size is increasedEl3 and their quantum number turn out to b e consistent 
with the constraints derived from the general RVB picture 
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6.4.5. Deconfined spinous 

The SL phase described above is expected to have deconfined spinons 
{S = \ excitations). These excitations should show up as an incoherent 
continuum in the spin-spin dynamical structure factor. However such a fea- 
ture would probably be rather difficult to observe on small 2D lattices, in 
particular due to the small number of inequivalent k-vectors in the Bril- 
louin zone. On the other hand, the binding energy of two spinons can be 
evaluated by comparing the ground-state energy and the first magnetic ex- 
citation energy on even and odd samples. In the case of the MSE model at 
J4 = 1, J2 = — 2 the results show the existence of a bound-state (it is more 
favorable to put two spinons in the same small sample than in separate 
ones, which is not surprising) but this does natmean that the spinons are 
confined (contrary to the conclusions of Ref. IHj). Interestingly this bind- 
ing energy seems to go to zero for the largest available sizes (Fig. 11): this 
might indicate the absence of attraction between spinons for large enough 
separation and an asymptotic deconfinement. 

It is important to stress here that the RVB SL phase discussed here 



(and its QDM counterparts of §5.5 and §5.6) is not the only way to spinon 
deconfinement in 2D. There is at least another scenario, inherited from one 
dimension, which is the sliding Luttinger liquid. Indeed, the Luttinger liquid 
behavior and the one-dimensional critical behavior of magnetic chains seem 
to be robust to small (or moderate) frustrating transverse couplings between 
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Fig. 11. Spinous biudiug euergy as a function of the system size in the MSE model 
at J2 = — 2 and J4 = 1. The vertical bars correspond to the range of values found for 
different sample shapes. 



chains, as o bser ved both in tbenrptir-a r^°l ' P^iH^^4 P^4 P^"^fP^5|P^^ and numerical 
approaches This regime bet ween one and two dimensions which may 

is the subject of a number of recent 



have beeiL XiLserved in CS2CUCI4 
studies.c2l'E2i 



6.5. Other models with MSE interactions 

Multiple-spin interactions are present in a number of models that were 
found to exhibit fractionalization or an RVB hquid ground-state. Well 
known examples of MSE interactions with an Ising symmetry are Z2 gauge 
theories, where the gauge invariant plaquette term is a product of P auli 
matrices Yii • Such theories have a deconfined phase in 2+1 dimcnsior££3 
and their relevance to fractionalized phases of 2D electronic systems has 



been pointed out by Senthil and Fisher.E£il The connexion between Z2 



gauge theories and QDM was mentioned in §5.6.7. Some MSE spin models 



with an I sing symmetry and a fr actio nalized ground-state were discussed 



by KitaevJ^Ha Nayak and Shtengel.E3 In the other Hmit of a f7(l) (or XY) 



symmetry several models have been studied. Recent e xamples are based 
upon the spin-i four-spin XY ring exchange interaction E§ ^ ' eI'^^ 



n = -K 



E 

{ijki) 



{s+s. 



H.c 



(49) 



which is the XY analog of the SU{2) MSE interaction P1234 + H.c. 
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7. Antiferromagnets on the kagome (and related) lattices 

In RVB SL phases described in j |5.5| , § |5.6| and § |6.4| singlet excitations ( vi- 
sons) are gapped. In the Z2 ga uge the ory approach this gap is essential 
for the consistency of the theoryJ^^'Enl The gauge-field quasi particles are 
vortices of the gauge field and carry a unit Z2 gauge flux but no spin. As 



illustrated in §5.6.6 in the framework of a simple QDM, visons have long- 
ranged interaction with spinous. If the spectrum of these visons has a gap 
then the spinous are unconfined and the phase is "fractionalized" . If they 
condense, the long range interaction between them and the spinous frus- 
trate the motions of the latter which remained confined. The gap in the 
singlet sector (above the topological degeneracy) is thus a crucial ingre- 
dient of these RVB SL. At the end of this section we will show how the 
spin-i first-neighbor Heisenberg model on the kagome lattice represents 
an enigma with respect to this scheme (as well as the QDM described in 



Sec. 5.7). We begin the section by a review of the properties of various 
models on the kagome lattice, which topology is the source of an extreme 
frustration, and, in many models, source of an extensive degeneracy of the 
ground-state. 

7.1. Miscellaneous models on the kagome lattice 

There has been a large number of studies devoted to different antiferro- 
magnetic models on the kagome lattice. The next neighbor Ising model on 
such a lattice is disordered, its entropy per site is very large S^^^ = 0.502, 
more than half the independent spin value, much larger than the triangu- 
lar lattice value S't™* = 0.323 and of the orde r of Pauling approximation 



for independent triangles ^Pauiing = 0.501.c23 This suggests that the cor- 



relations in this systern are very weak: the model remains disordered at 



all temperatures .c23c£3 Moessner and Sondhi have studied this Ising model 



in a transverse magnetic field (the simplest way to include some quantum 
fluctuations in the rnorjel) the model fails to order for any transverse field, 
at any temperature. 



145 



237 



The nearest neighbor classical Heisenberg model on the kagome lattice 
has also a huge ground-state degeneracy. This property is easy to under- 
stand and holds on different lattices with corner sharing units such as the 
checkerboard l attic e or the three dimensional pyrochlore lattice (Moessner 
and Chalkerli£3't23). On all these lattices the nearest neighbor Heisenberg 
Hamiltonian can be written as the sum of the square of the total spin Sa of 
individual units a (a tetrahedron in the 2-d and 3-d pyrochlore cases and 
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a triangle for the kagome lattice), which share only one vertex. Classical 
ground-states are obtained whenever Va Sa = 0. This condition fixes the 
relative positions of the three classical spins of a triangle at 120 degrees 
from each other in a plane. But it does not fix the relative orientation of 
the plane of a triad with respect to the planes of triads on corner sharing 
triangles: the model has a continuous local degeneracy ^ ^ at T = 0." 
Thermal fluctuations select coplanar configurations E23'E£3'E13 The possibil- 
ity of long-range order in spin-spin correlations at vrn-y lo w temperature has 
been discussed without any definitive conclusion.E23'EiJ The order param- 
eter of the planar phase is defined by the local helicity (sometimes called 
vectorial chirality) : 



C = 5i X 52 + 52 X ^3 + ^3 X 5i 



(50) 



where the three sites define a triangle. This kind of order is sometimes called 
nematic by analogy to liquid crystals. The existence of such an order param- 
eter might be related to the instability of the classical Heisenberg model on 
the kagome lattice to Dzyaloshinsky-Moriya interactions .tlj'Efj The clagai 
cal model has a large density of low- lying excitations at low temperature 



7.2. Spin-^ Heisenberg model on the kagome lattice: an 
extreme play-ground for "quantum fluctuations" 

The nearest neighbor spin-i quantum Heisenberg model on the kagome 

object of many 



lattice 
studiesElEr 



From these studies the follow- 



ing important facts have emerged : 



7.2.1. Ground-state energy per spin 

The Heisenberg model on the kagome lattice has an extremely low energy 
per bond ((S^.S^) ~ —0.44) ^ 87% of the energy per bond in an isolated 
triangle. On this lattice the energy per bond of the spin-i system is much 
lower than the classical energy ^ 1.74, a ratio much larger than in any 
other 2D magnet, that can only be compared to the value obtained for the 
Bethe chain (1.77). The kagome lattice is the 2D lattice which offers the 
largest stabilization due to quantum fluctuations. 



"Counting the planar ground-states amounts to determine in how many ways one can 
associate one of the three letters A, B and C to each site so that each triangle ha p iipina 
along the three different orientations. This already represents an extensive entropyC^jL^J 
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7.2.2. Correlations 



The ground-state is disorderedE^ Within acmiracy of the finite-size nu- 



merical computations, spin-spin correl; 
(Fig. n2[), chirahty-chirahty correlations 



iio 



246 



ns,cl3 dimer-dimer correlations 
are short-rangecLwhich is con- 



sistent with the previous point and series expansion results.!^ 
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Fig. 12. Dimer-dimer correlations in the ground-sta te o f the kagome Heisenberg model 
(black triangles) and in the MSE model discussed in § 3.4 (open symbols) versus distance. 
Although the decrease of these correlations is weaker in the KH model than in the MSE 
model, it seems p#"*ef*heless roughly exponential in the first two decades, as the spin-spin 
correlations are.Ef3E13 



7.2.3. Spin- gap 

There is nlausibly a spin-gap of the order of 1/20*'' of the coupling 
constant.E£3 In view of the smallness of this spin-gap with regards to the 
available sizes caution is necessary. The above conclusion was drawn from 
the value of the microscopic spin gap Eq(S =^af in + 1) ~ Eo{S — Smin) on 
exact spectra of samples with up to 36 spins,E£2l where Eq{S) is the lowest 
energy in the S sector, and Smin = or i depending on the parity of the 
total number of spins. The finite-size effects on these results are an order 
of magnitude smaller than in a Neel ordered antiferromagnet. Nevertheless 
they are still not negligible for these sizes. An alternative determination 
of the spin-gap can be obtained along the following line. The lowest exact 
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eigenstate in each total spin sector S" of a spins sample defines the en- 
ergy per spin e at T = as a function of its magnetization m = S/{N/2). 
For low positive values of the magnetization, one can fit e(m) to the phe- 
nomenological law:^ 

e(m) = e(0) + am + 6mV2 + 0(m3) (51) 

a and b are intensive quantities which depend on but should con- 
verge to some (possibly zero) value when iV — > oo. This point of view 
is important because for an infinite system the thermodynamic function 
e(m) = ImiN^ac E(S — mN/2)/N can be measured (through the zero 
temperature magnetization curve), but not E{S 0(1)). The physical 
significance of Eq. ^ is clear: a measures half the spin gap (in a thermody- 
namic sense, that is the magnetic field required to magnetize the system 
at zero temperature) and 

where x is the homogeneous susceptibility of the medium for fields larger 
than the critical field He — a. This determination jiLthe thermodynamic 
spin gap leads to a renormalization of the raw datac^ for small sizes (see 
Fig. ^3|). A linear extrapolation versus 1/N (which should give a lower 
bound of the spin-gap) leads to the value 0.06 for the spin gap (coo = 
—0.4365 and Xoo — 0.34). This determination is in agreement with the direct 
extrapolation of the microscopic gap Eo{S = Smin + 1) — Eq[S = Smin)- 
For = 36 at the smallest non-zero magnetization, the linear term of 
equation (^ is 90% of the quadratic term: this is an estimate of the degree 
of confidence on the existence of a spin-gap. 

7.2.4. An exceptional density of low lying excitations in the singlet 
sector 

Whatever the ultimate fate of the spin gap a big surprise emerges from 
the exact spectra: the probable absence of gap in the singlet sector and the 
anomalous density of low energy states adjacent to the ground-state. Let us 
first comment the second point in details: even on the smallest size spectra 
the low lying states appear contiguous to the ground-state and the spectra 



^This phenomenological form cannot extend below m = and beyond m = 1/3 because 
at both magnetizations an angular point app nara i in l e (m) with a discontinuity of the first 
derivative signalinp-p. magnetization plateau. t£ZrE£3 This is discussed in more details in 
our lectures notes. 
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Fig. 13. Finite size scaling of the spin gap in the spin-^ Heisenberg model "H = 
S<i j> Si.Sj on the kagome lattice. Triangles (rosp. squares) are the raw results of 
the microscopic spin gap in the even (rosp. odd) samples. Bullets represent the "ther- 
modynamic" spin gap, they are obtained by the procedure described in the text (errors 
bars come from the rms uncertainty in the fits). 



are extremely dense. The number of singlet levels in the spin-gap (taken 

as a natural energy band-width of the problem) increases exponentially 

fast with iV as ~ 1.15^ as far as the < 36 systems are concerned (see 

Fig. This property remains (and can be checked on larger systems) 

when the Hilbert space is restricted to that of first-neighbor valence-bond 
25| i7q 



coverings 




Fig. 14. Logarithm of the number of singlet states in the spin gap versus sample size 
(black squares). The short dashed and long dashed curves display the theoretical law 
(Eq. |54| with d = 2) (short dashes: p = 1, long dashes: p = 2). 
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Some remarks are necessary to fully appreciate this property. The 2^ 
states of the system are stretched on an energy scale of the order of NJ 
where J is the coupling constant of the Hamiltonian. This implies that on 
most of the spectrum the density of states increases exponentially with N. 
If we specialize to the 5 = sector as we will do below, the picture is not 
very different: the number of states is — ~ ^(iv") ^'^d here too, 

in most of the spectrum the density is exponentially increasing with N . 
But in all the phases that we have studied up to now, the nature of the 
ground-state and of the low lying excitations leads to a different behavior 
at the bottom of the spectrum. The ground-state degeneracy is 0(1) in 



VBC, VBS and in the RVB SL (of the type discussed in |5j and §|6j for 
instance) and it is 0{NP) in Neel ordered states with p sublattices. In all 
these situations the low-lying excitations are described as modes or quasi- 
particles. Let us suppose that single-particle excitations have a dispersion 
law: 

a;(k) - fcP. (53) 

In a d-dimensional system the total internal energy increases with tem- 
perature as E ^ NT'^p+d)/p^ the specific heat as Cy - NT'^/p as well as 
the entropy S - NT'^^p. In a micro-canonical point of view the density of 
states p{E) or a large enough system is simply related to the entropy by 
S = \og{p{E)). Since S - NT'^^p = N{E/NY^^p+'^'> we obtain 

\og{p{E)) cx N-!^ (54) 

As an example, let us consider the Rokhsar Kivelson QDM on the square 
lattice at the quantum critical point (Sec. ||): the dispersion law of the so- 
called resonons is quadratic around Q = (tt, tt), the logarithm of the number 
of states increases as (long dashes of Fig. |l|). 

Even with such many-particle excitations one would expect a number of 
levels increasing more slowly than in the numerical spectra. Infinitely soft 
low-energy modes (p ^ cx)), are necessary to recover a density of low- lying 
levels growing as ~ . It is still unclear if we can do a connection between 
the "zero modes" of the classical model at T = and this picture. And we 
cannot completely indulge ourselves in saying that quantum fluctuations 
are unable to lift the classical degeneracy as quantum fluctuations seem to 
open a spin gap. 

A physical consequence of this exceptional density of low lying sin- 
glets can be observed in the specific heat: at low temperature the spe- 
cific heat of this spin system is unusually large, with a double peak 



February 2, 2008 4:30 



WSPC/Trim Size: 9in x 6in for Review Volume 



a 



64 



G. Misguich and C. Lhuillier 



252 



structure,^' El ?£^'^ insensitive to relatively large magnetic fields 
This is easily understood if we sunpose that in this energy range there 
is a large density of singlet states.c^ T,lii^ result is to be compared to 



the experimental results of Ramirez et aZ££3 on SrCrgpGai2-9pOi9 (notice 



however that it is a spin-| compound and that numerical calculations were 
performed on the spin-i model) where the specific heat around 5 K has an 
extremely low sensitivity to magnetic fields up to 10 Tesla, whereas the 
homogeneous susceptibility in this range of tenmCTature is probably very 
low if we notice that it turns down around 50-ftrH£3 



7.2.5. Absence of gap in the singlet sector 



Up to = 36 sites there is no evidence of a possible gap in the singlet 
spectrum: this is an exceptional phenomenon in quantum mechanics of 



y the rule 
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that the ground-state of 



small systems where discretization is usual!; 

It has been advocated in recent papcrsL£ff 
this model could break the translational symmetry and be a VBC. The first 
proposed crystalE£j is made of resonating stars with 6 dimers. The corre- 
sponding unit cell has 12 spins. The second VBCEHj is made of resonating 
(trimerized) hexagons and was already discussed as the most reasonable 
crystal (if any is realized in this system, whi ch w e think is not clear at all) 



by Marston and ZengElJ and Zeng and Elser.EIil This later VBC has a unit 



cell with 36 sites. In both scenarios it is the energy gain obtained by lo- 
cal resonances (involving respectively 6 and 3 valence-bonds) which drives 
the system toward a VBC. From the energy point of view the star VBC is 
(by far) less realistic since it involves a much longer resonance loop.^ This 
resonance loop involving 6 valence-bonds around a star has a vanishing am- 
plitude at the lowest non-trivial order of the overlap expansion in the RVB 
subspace, as was shown by Zeng and Elseillld. To our opinion there is no 
reason to think that the physics of the kagome model can be described from 
the limit of weakly coupled stars. On the other hand, in the approximation 
where only the shortest resonance loops are present, the model was indeed 
found to be in the hexagon VBC phase, as re-discovered recently frorna 
different point of viewH£3 A crucial (numerical) result of Zeng and Elsenill 
is however that this VBC melts when higher order resonances loops are 
included. 



"'Notice aiseithat the gap associated to the star VBC is predicted to be of the order of 
1/10 of Jc£j which is 10 times larger than the largest distance between two consecutive 
levels in the spectrum of the Af = 36 sample. 
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We have studied the 36-sites sample which can accommodate these two 
VBC: the low-lying levels of the spectrum do not give a clear picture of 
the supposed to be VBC. The eigenstates with the quantum numbers cor- 
responding to these two VBC are not the lowest energy levels in the spec- 
trum. No gap is seen in the spectrum, the largest distance between two 
consecutive states is 10~^ and seems distributed at random, whereas the 
average distance between two consecutive states in the 50 first states is: 
2.510~^. It has been argued that long wave-length quantum fluctuations 
(almost absent in the A'^ = 36 sample which contains only two resonat- 
ing hexagons or three stars) could eventually restore the order. We think 
that this is incorrect. In a VBC it is local resonances which favor the crys- 
tal and longer resonances which tend to reduce the order parameter and 
which could eventually destabilize it. Long- wavelength will tend to restore 
a larger ground-state symmetry (reducing the degeneracy). From this point 
of view it is unlikely that a VBC pattern will appear in larger systems if 
it is not apparent in the smallest systems (provided that boundary condi- 
tions do not frustrate the corresponding VBC). We thus consider that up 
to now numerical results do not support the claims of a star or hexagon 
VBC ground-state in this model. 

If the hexagon VBC was however realized the important question would 
probably be why the associated gap is so tiny, certainly smaller than 10^"^. 
A natural scenario would be that of a very close proximity to a quantum 
critical point, but that critical point remains to be identified. A third pos- 
sibility would be that of a system with a true extensive entropy at zero 
temperature. This is less likely in our opinion since (up to now) only mod- 
els with simple (if not trivial) local symmetries were found to have such 
ground-state entropy (Ising antiferromagnet on the triangular lattice for 



instance or the QDM of §5.7) and this entropy is usually sent to (possibly 



small but) finite temperatures by almost any perturbation. 



7.2.6. Anomalous density of states in other spin sectors 

An anomalous density of low lying states is equally observed in the spin ^ 
sector (where the law could be fitted to A^l.15^), in the spin 1 sector as well 
as in other sectors with larger total spin. It should be noticed that such a 
density of states implies the absence of an intrinsic energy scale for the low 
lying excita,tions a, phenomenon that has been observed in inelastic neutron 



scattering (Ref. lH3 and Refs. therein) and theoretically in the imaginary 



part of the dynamic susceptibility calculated within the dynamical mean 
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field theory (Georges et a/.E£ZI). A hig±i_spin susceptibility just above the 



266 



spin gap is not excluded in the spin- 2 ' 

The global picture of this phase is thus that of a SL with no long- 
ranged correlations in any local observable, and an large entropy of singlets 
at T <C J, which is a manifestation of the extraordinary large density of 
states in each S subspace. 



7.3. Next-neighbor RVB description of the spin-^ kagome 
antiferromagnet 



Considering a supposed-to-be large spin-gap, Zeng and Elscnilil proposed 
a description of the ground-state and low-lying excitations of the kagome 
model in the basis of next neighbor valence bonds. They analyzed in this 
context the dimer dynamics and showed on a = 36 sample that the 
hexagon VBC -favored by the shortest (three-)dimer m oves - m elts when 



introducing higher order tunneling. Mila and Mambrinic£il'E£3 confirmed 



that this reduced Hilbert space of next neighbor valence bonds captures 
some of the most perplexing features of this magnet and specifically the 
absence of (measurable) gap in the singlet sector and the exponential num- 
ber of singlets. One of us, D. Serban and V. Pasquiei£Z3 have elaborated 
on this work and on Zcng and Elser's approach and proposed a QDM with 
an extensive zero-point entropy and critical (energy-energy'^) correlations 
(see §5/7). All these results point to an absence of an intrinsic low-energy 
scale. This feature is typical of a critical state, but as seen in the above 
discussion, the simple RK picture does not seem to fit nicely to the exact 
diagonalization data: may be the available sizes are too small or the be- 
havior of this quantum system_corresponds to something completely new. 
Some recent numerical result£^ (in the full spin-i Hilbert space as well as 
in the RVB subspace) showed that (static) non-magnetic impurities (holes) 
experience an unexpected repulsion in this system and that no significant 
magne 
argued 



in 
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moment is created in the vicinity of the impurities. It has been 
that static non-magnetic impurities are useful to detect a possi- 
ble spinon deconfinement in two-dimensional antiferromagnets. From this 
point of view the results mentioned above suggest such a deconfinement. 



^Dimer-dimer correlations are short-ranged in this model. 
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7.4. Experiments in compounds with kagome-like lattices 

The low temperature specific heat of SrCrgGai20i9 (the magnetic chromium 



ion has a spin 3/2) is apparently dominated by singlet statesE£3 The mag 



netic excitations of this compound as_aeen by muons can be described as 



spins ^ itinerant in a sea of singlets. clill The non-linear spin susceptibility 



of SrCr9Gai20ifl_exhibits a very large increase at about 5 K, reminiscent 
of spin glasses,EId but neutrons and muons show that a very significant 
fraction of the spins are n ot fr ozen below this temperature and exhibit 
still very rapid fluctuations. EI3 The same phenomena have been observed 
in two jarosites that are equally good rno dels of kagome antiferromagnets 
with half-odd-integer spin per unit cell 



i 



273,274 



7.5. "Haldane's conjecture" 

Whereas the classical Heisenberg model on the kagome, checkerboard and 
pyrochlore lattices share the properties of local continuous degeneracy and 
disorder at T = 0, their quantum counterparts are quite different. As it has 
been explained in Sec. ^, the Heisenberg model on the checkerboard lattice 
has an ordered VBC with gaps to all excitations. In contrast to the case 
of the 5' = i Heisenberg model on the kagome lattice, recent results from 
Hida£3 show that there is a large gap to all excitations in the 5* = 1 Heisen 
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berg model on this same lattice, in agreement with experiments 

Less is known on the ground-state of the Heisentecg model on the three- 
dimensional pyrochlore lattice: Canals and LacroixEZa have shown that the 
spin-spin correlations are short ranged and they have observed that on a 
16-sites spectrum the first excitations are singlet ones. _Eouet et al. have 
computed the spectrum of a 32-sites pyrochlore sample c3, this work con- 
firms that the first excitations are still singlets for this size. There is plau- 
sibly a degeneracy of the ground-state in the thermodynamic limit but 
no evidence of a closing of the gap above: the system may therefore be 
a VBC. A first descrmlion of the singlet sector was proposed by Harrip . 
Berlinsky and BrudeiEZa and was developed further by Tsunetsugri223^E£l . 



It starts from the limit of weakly-coupled tetrahedron (and thus breaks 
some spatial symmetries). The ground-state of an isolated tetrahedron is 
two-fold degenerate and can be described by an Ising pseudo-spin. An ef- 
fective Ha,r niltp nia,n describing the interactions between these pseudo-spins 
was writtenEIj and analyzed in a semi-classical approximation. It was con- 
cluded that the system maybe a VBC. It wasm-gued by Tsunetsugu that 



a soft mode could exist in the singlet sector.Eill. The CORE approach of 
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Berg et aZjliB seems more appropriate to deal with these systems where 
dimerization is probably the dominant phenomenon. They treated larger 
units (block of four tetrahedron) and concludes to t he pyis tence of a VBC 
with a larger unit cell than the one predicted befora2I3l2^'E^ and a small 
singlet gap. 

All these results seem to confirm a 2D version of Haldane's conjecture: 
among these frustrated systems with local continuous degeneracies in the 
classical limit, the spin-^ kagome antiferromagnet is the only system to 
have an half-odd integer spin in the unit cell. It is maybe not by chance 
that it is the only one with gapless excitations.^ The spin-i Heisenberg 
model on the checkerboard lattice or on the pyrochlore lattice and the spin- 
1 Heisenberg model on the kagome lattice have integer spins in the unit cell 
and quantum fluctuations lead to gapful excitations. It is also interesting 
to note that these results pr q co n si s tent with a generalization of the LSM 



216,217 



157 2 



£n 

21g| 



21£ 



theorem in dimension two 

An interesting analysis of Von Delft and IIenlejc£3 supports this con- 
jecture. These authors studied the collective tunneling of a small cluster 
of spins between two spin configurations that are degenerate in the clas- 
sical limit. They found that for half-odd-integer spins the tunneling am- 
plitude for a cluster of six spins around an hexagon (and on other larger 
loops) is exactly zero because of destructive interferences between different 
symmetry-related instantons. For small integer spins the interference is con- 
structive and the tunnel amplitude and the tunnel splitting are large: this 
is consistent with numerical results which gives a large gap for the 5 = 1 
kagome antiferromagnet, and small gaps (if any) in the spin-^ system. 



8. Conclusions 

We conclude by summarizing some properties -and related open questions- 
of the different phases discussed in this review. 

The properties of these phases are summarized in Table ^. Semi-classical 
phases with Neel long ranged order, magnons as gapless excitations, do exist 
in spin-i 2D systems with moderate frustration: the Heisenberg model on 
the triangular lattice is the most explicit example, with a sublattice magne- 
tization about one half of the classical value.i The ground-state degeneracy 
is in the thermodynamic limit a power of N . 



''The SCGO compound wi|th-7 spins 3/2 by unit cell belongs to the same category and 
as also gapless excitations. t£j 
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An increased frustration, lower coordination number or smaller spin lead 
to quantum phases, with a ground-state of higher symmetry, no long ranged 
order in spin-spin correlations, a spin gap and the restored SU (2) symmetry. 
Two main alternatives are then opened: the VBC or VBS phases on one 
hand, the RVB SL on the other. RVB SL and VBC (as well as VBS) first 
requires the formation of local singlets. When a particular local resonance 
pattern dominates the dynamics of the Hamiltonian the system will try to 
maximize the number of occurrence of this pattern. This is usually achieved 
by a regular arrangement, that is a VBC. When no such pattern dominates 
the system may form a translation invariant RVB SL. In the first case the 
ground-state can be qualitatively described by one ordered configuration 
of singlets dressed by small fluctuations. In the RVB SL the amplitudes 
of the wave-function are distributed over an exponentially large number of 
configurations. These ground-states lead to very different excitations: AS* — 
1 gapped magnons in the first case (and A^* — domain- wall excitations), 
gapped AS — visons and gapped AS* — ^ unconfincd spinous in the 
second case. 

In agreement with the large- iV results (Sec. ^), VBC or VBS phases 
appear in general in quantum situations where the large-S* classical limit 
displays coUinear order ^, whereas up to now RVB SL phases have only 
be encniintere d in ran ge of parameters were the classical solutions are non 
collinear (MSEllMlliJ and J1-J2 on the honeycomb latticeEZl). 

These states obey the 2D extension of LSM theorem: if 25' is odd in the 
unit cell and if excitations are gapped there is a ground-state degeneracy in 
the thermodynamic limit (with periodic boundary conditions) . However the 
origin of the degeneracy differs in the two types of quantum phases. In the 
VBC phases the degeneracy is associated to spontaneously broken trans- 
lation symmetry whereas in the RVB SL the degeneracy has a topological 
origin. In the VBS (or explicit VBC) the ground-state is unique but 25 is 
even in the unit cell. As in one dimension the LSM theorem seems to play 
an important role and a, formal proof of its validity in 2D is perhaps very 



close, if not achieved.EJ A RVB SL may have been observed numerically 
in a spin-^ model on the hexagonal latticcEH From the point of view of 
large- and QDM approaches a topological degeneracy is expected in such 
a SL. Interestingly that degeneracy is not imposed by the 2D extension of 
LSM's theorem and has not been detected.0 

These paradigms are relatively well understood, at least on the qualita- 



'^We have discussed some possible counter examples in section 3. 6 
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tive level. They also appjear naturally in the broader context of the classifica- 
tion of Mott insulators .c3 However several kinds of 2D frustrated magnets 
do not fall in these simple classes and many open questions remain. 

This review was restricted to SU{2) invariant Hamiltonians. Whereas 
the Ising limit has been much studied, the differences between quantum 
XY and Heiscnberg models have received much less attention. 

Chiral SL have not been discussed in this review. They are character- 
ized by a broken time- reversal syTTiiTietry. This possibility has been studied 
intensely since the 80'sEME£j'E£j'Q. To the best of our knowledge it has not 
yet been identified in a realistic 2D model.'' 

The issue of quantum phase transitions in frustrated antiferromagnet is 
also an active topic that is not presented in this review. Many properties 
of these critical points are still unknown, not to mention the fascinating 
problems associated with (quenched) disorder. 

Limited by place (and competence) we have not discussed in details the 
works done on spatially anisotropic models. This field which is in between 
one (review by P. Lecheminant in this book) and 2D is extremely fiourishing 
tackled by bosonization and large-iV methods. 

To conclude we would like to emphasize that new analytical and/or 
numerical methods are highly desirable to proceed in the analysis of the 
two emblematic problems by which we have opened and closed this review: 
the J1-J2 model on the square lattice and the spin-i Heisenberg model on 
the kagome (and pyrochlore) lattices. In both of these problems a consensus 
remains to be obtained. 
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